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Abstract. We use numerical climate simulations, paleocli- 
mate data, and modern observations to study the effect of 
growing ice melt from Antarctica and Greenland. Meltwa- 
ter tends to stabilize the ocean column, inducing amplifying 
feedbacks that increase subsurface ocean warming and ice 
shelf melting. Cold meltwater and induced dynamical effects 
cause ocean surface cooling in the Southern Ocean and North 
Atlantic, thus increasing Earth’s energy imbalance and heat 
flux into most of the global ocean’s surface. Southern Ocean 
surface cooling, while lower latitudes are warming, increases 
precipitation on the Southern Ocean, increasing ocean strat- 
ification, slowing deepwater formation, and increasing ice 
sheet mass loss. These feedbacks make ice sheets in contact 
with the ocean vulnerable to accelerating disintegration. We 


hypothesize that ice mass loss from the most vulnerable ice, 
sufficient to raise sea level several meters, is better approx- 
imated as exponential than by a more linear response. Dou- 
bling times of 10, 20 or 40 years yield multi-meter sea level 
rise in about 50, 100 or 200 years. Recent ice melt doubling 
times are near the lower end of the 10^-0-year range, but 
the record is too short to confirm the nature of the response. 
The feedbacks, including subsurface ocean warming, help 
explain paleoclimate data and point to a dominant Southern 
Ocean role in controlling atmospheric CO 2 , which in turn ex- 
ercised tight control on global temperature and sea level. The 
millennial (500-2000-year) timescale of deep-ocean ventila- 
tion affects the timescale for natural CO 2 change and thus 
the timescale for paleo-global climate, ice sheet, and sea 
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level changes, but this paleo-millennial timescale should not 
be misinterpreted as the timescale for ice sheet response to 
a rapid, large, human-made climate forcing. These climate 
feedbacks aid interpretation of events late in the prior in- 
terglacial, when sea level rose to +6-9 m with evidence of 
extreme storms while Earth was less than 1 °C warmer than 
today. Ice melt cooling of the North Atlantic and Southern 
oceans increases atmospheric temperature gradients, eddy 
kinetic energy and baroclinicity, thus driving more power- 
ful storms. The modeling, paleoclimate evidence, and on- 
going observations together imply that 2 °C global warming 
above the preindustrial level could be dangerous. Continued 
high fossil fuel emissions this century are predicted to yield 
(1) cooling of the Southern Ocean, especially in the Western 
Hemisphere; (2) slowing of the Southern Ocean overturning 
circulation, warming of the ice shelves, and growing ice sheet 
mass loss; (3) slowdown and eventual shutdown of the At- 
lantic overturning circulation with cooling of the North At- 
lantic region; (4) increasingly powerful storms; and (5) non- 
linearly growing sea level rise, reaching several meters over 
a timescale of 50-150 years. These predictions, especially 
the cooling in the Southern Ocean and North Atlantic with 
markedly reduced warming or even cooling in Europe, dif- 
fer fundamentally from existing climate change assessments. 
We discuss observations and modeling studies needed to re- 
fute or clarify these assertions. 


1 Introduction 

Humanity is rapidly extracting and burning fossil fuels with- 
out full understanding of the consequences. Current assess- 
ments place emphasis on practical effects such as increas- 
ing extremes of heat waves, droughts, heavy rainfall, floods, 
and encroaching seas (IPCC, 2014; USNCA, 2014). These 
assessments and our recent study (Hansen et ah, 2013a) con- 
clude that there is an urgency to slow carbon dioxide (CO 2 ) 
emissions, because the longevity of the carbon in the climate 
system (Archer, 2005) and persistence of the induced warm- 
ing (Solomon et ah, 2010) may lock in unavoidable, highly 
undesirable consequences. 

Despite these warnings, fossil fuels remain the world’s pri- 
mary energy source and global CO 2 emissions continue at 
a high level, perhaps with an expectation that humanity can 
adapt to climate change and find ways to minimize effects via 
advanced technologies. We suggest that this viewpoint fails 
to appreciate the nature of the threat posed by ice sheet insta- 
bility and sea level rise. If the ocean continues to accumulate 
heat and increase melting of marine-terminating ice shelves 
of Antarctica and Greenland, a point will be reached at which 
it is impossible to avoid large-scale ice sheet disintegration 
with sea level rise of at least several meters. The economic 
and social cost of losing functionality of all coastal cities is 
practically incalculable. We suggest that a strategy relying 


J. Hansen et al.: Ice melt, sea level rise and superstorms 

on adaptation to such consequences will be unacceptable to 
most of humanity, so it is important to understand this threat 
as soon as possible. 

We investigate the climate threat using a combination of 
atmosphere-ocean modeling, information from paleoclimate 
data, and observations of ongoing climate change. Each of 
these has limitations: modeling is an imperfect representa- 
tion of the climate system, paleo-data consist mainly of proxy 
climate information usually with substantial ambiguities, and 
modern observations are limited in scope and accuracy. How- 
ever, with the help of a large body of research by the scientific 
community, it is possible to draw meaningful conclusions. 

2 Background information and organization of the 
paper 

Our study germinated a decade ago. Hansen (2005, 2007) ar- 
gued that the modest 21st century sea level rise projected by 
IPCC (2001), less than a meter, was inconsistent with pre- 
sumed climate forcings, which were larger than paleoclimate 
forcings associated with sea level rise of many meters. His 
argument about the potential rate of sea level rise was nec- 
essarily heuristic, because ice sheet models are at an early 
stage of development, depending sensitively on many pro- 
cesses that are poorly understood. This uncertainty is illus- 
trated by Pollard et al. (2015), who found that addition of 
hydro-fracturing and cliff failure into their ice sheet model 
increased simulated sea level rise from 2 to 17 m, in response 
to only 2 °C ocean warming and accelerated the time for sub- 
stantial change from several centuries to several decades. 

The focus for our paper developed in 2007, when the 
first author (JH) read several papers by co-author P. Hearty. 
Hearty used geologic field data to make a persuasive case for 
rapid sea level rise late in the prior interglacial period to a 
height +6-9 m relative to today, and he presented evidence 
of strong storms in the Bahamas and Bermuda at that time. 
Hearty’s data suggested violent climate behavior on a planet 
only slightly warmer than today. 

Our study was designed to shed light on, or at least raise 
questions about, physical processes that could help account 
for the paleoclimate data and have relevance to ongoing and 
future climate change. Our assumption was that extraction 
of significant information on these processes would require 
use of and analysis of (1) climate modeling, (2) paleoclimate 
data, and (3) modern observations. It is the combination of all 
of these that helps us interpret the intricate paleoclimate data 
and extract implications about future sea level and storms. 

Our approach is to postulate existence of feedbacks that 
can rapidly accelerate ice melt, impose such rapidly grow- 
ing freshwater injection on a climate model, and look for 
a climate response that supports such acceleration. Our im- 
posed ice melt grows nonlinearly in time, specifically expo- 
nentially, so the rate is characterized by a doubling time. To- 
tal amounts of freshwater injection are chosen in the range 
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1-5 m of sea level, amounts that can be provided by vulner- 
able ice masses in contact with the ocean. We find signif- 
icant impact of meltwater on global climate and feedbacks 
that support ice melt acceleration. We obtain this informa- 
tion without use of ice sheet models, which are still at an 
early stage of development, in contrast to global general cir- 
culation models that were developed over more than half a 
century and do a capable job of simulating atmosphere and 
ocean circulation. 

Our principal finding concerns the effect of meltwater on 
stratification of the high-latitude ocean and resulting ocean 
heat sequestration that leads to melting of ice shelves and 
catastrophic ice sheet collapse. Stratification contrasts with 
homogenization. Winter conditions on parts of the North At- 
lantic Ocean and around the edges of Antarctica normally 
produce cold, salty water that is dense enough to sink to the 
deep ocean, thus stirring and tending to homogenize the wa- 
ter column. Injection of fresh meltwater reduces the density 
of the upper ocean wind-stirred mixed layer, thus reducing 
the rate at which cold surface water sinks in winter at high 
latitudes. Vertical mixing normally brings warmer water to 
the surface, where heat is released to the atmosphere and 
space. Thus the increased stratification due to freshwater in- 
jection causes heat to be retained at ocean depth, where it 
is available to melt ice shelves. Despite improvements that 
we make in our ocean model, which allow Antarctic Bottom 
Water to be formed at proper locations, we suggest that ex- 
cessive mixing in many climate models, ours included, lim- 
its this stratification effect. Thus, human impact on ice sheets 
and sea level may be even more imminent than in our model, 
suggesting a need for confirmatory observations. 

Our paper published in Atmospheric Chemistry and 
Physics Discussion was organized in the chronological order 
of our investigation. Here we reorganize the work to make 
the science easier to follow. First, we describe our climate 
simulations with specified growing freshwater sources in the 
North Atlantic and Southern oceans. Second, we analyze pa- 
leoclimate data for evidence of these processes and possible 
implications for the future. Third, we examine modern data 
for evidence that the simulated climate changes are already 
occurring. 

We use paleoclimate data to find support for and deeper 
understanding of these processes, focusing especially on 
events in the last interglacial period warmer than today, 
called Marine Isotope Stage (MIS) 5e in studies of ocean 
sediment cores, Eemian in European climate studies, and 
sometimes Sangamonian in US literature (see Sect. 4.2 for 
timescale diagram of marine isotope stages). Accurately 
known changes of Earth’s astronomical configuration altered 
the seasonal and geographical distribution of incoming ra- 
diation during the Eemian. Resulting global warming was 
due to feedbacks that amplified the orbital forcing. While 
the Eemian is not an analog of future warming, it is useful 
for investigating climate feedbacks, including the interplay 
between ice melt at high latitudes and ocean circulation. 


3 Simulations of 1850-2300 climate change 

We make simulations for 1850-2300 with radiative forcings 
that were used in CMIP (Climate Model Intercomparison 
Project) simulations reported by IPCC (2007, 2013). This al- 
lows comparison of our present simulations with prior stud- 
ies. First, for the sake of later raising and discussing funda- 
mental questions about ocean mixing and climate response 
time, we define climate forcings and the relation of forcings 
to Earth’s energy imbalance and global temperature. 

3.1 Climate forcing, Earth’s energy imbalance, and 
climate response function 

A climate forcing is an imposed perturbation of Earth’s en- 
ergy balance, such as change in solar irradiance or a ra- 
diatively effective constituent of the atmosphere or surface. 
Non-radiative climate forcings are possible, e.g., change in 
Earth’s surface roughness or rotation rate, but these are small 
and radiative feedbacks likely dominate global climate re- 
sponse even in such cases. The net forcing driving climate 
change in our simulations (Fig. S16 in the Supplement) is 
almost 2 W m -2 at present and increases to 5-6 W m -2 at 
the end of this century, depending on how much the (nega- 
tive) aerosol forcing is assumed to reduce the greenhouse gas 
(GHG) forcing. The GHG forcing is based on IPCC scenario 
A1B. “Orbital” forcings, i.e., changes in the seasonal and ge- 
ographical distribution of insolation on millennial timescales 
caused by changes of Earth’s orbit and spin axis tilt, are 
near zero on global average, but they spur “slow feedbacks” 
of several W m , mainly change in surface reflectivity and 
GHGs. 

When a climate forcing changes, say solar irradiance in- 
creases or atmospheric CO 2 increases. Earth is temporarily 
out of energy balance, that is, more energy coming in than 
going out in these cases, so Earth’s temperature will increase 
until energy balance is restored. Earth’s energy imbalance is 
a result of the climate system’s inertia, i.e., the slowness of 
the surface temperature to respond to changing global cli- 
mate forcing. Earth’s energy imbalance is a function of ocean 
mixing, as well as climate forcing and climate sensitivity, the 
latter being the equilibrium global temperature response to a 
specified climate forcing. Earth’s present energy imbalance, 
+0.5-1 W m -2 (von Schuckmann et al., 2016), provides an 
indication of how much additional global warming is still 
“in the pipeline” if climate forcings remain unchanged. How- 
ever, climate change generated by today’s energy imbalance, 
especially the rate at which it occurs, is quite different than 
climate change in response to a new forcing of equal mag- 
nitude. Understanding this difference is relevant to issues 
raised in this paper. 

The different effect of old and new climate forcings is im- 
plicit in the shape of the climate response function, R(t), 
where R is the fraction of the equilibrium global temper- 
ature change achieved as a function of time following im- 


www.atmos-chem-phys.net/16/3761/2016/ 


Atmos. Chem. Phys., 16, 3761-3812, 2016 



3764 

position of a forcing. Global climate models find that a 
large fraction of the equilibrium response is obtained quickly, 
about half of the response occurring within several years, 
but the remainder is “recalcitrant” (Held et al., 2010), requir- 
ing many decades or even centuries for nearly complete re- 
sponse. Hansen (2008) showed that once a climate model’s 
response function R is known, based on simulations for an 
instant forcing, global temperature change, T it), in response 
to any climate forcing history, F(t), can be accurately ob- 
tained from a simple (Green’s function) integration of R over 
time: 

T(t) = j R(t)\dF/dt]dt. (1) 

dF / dt is the annual increment of the net forcing and the in- 
tegration begins before human-made climate forcing became 
substantial. 

We use these concepts in discussing evidence that most 
ocean models, ours included, are too diffusive. Such exces- 
sive mixing causes the Southern and North Atlantic oceans in 
the models to have an unrealistically slow response to surface 
meltwater injection. Implications include a more imminent 
threat of slowdowns of Antarctic Bottom Water and North 
Atlantic Deep Water formation than present models suggest, 
with regional and global climate impacts. 

3.2 Climate model 

Simulations are made with an improved version of a coarse- 
resolution model that allows long runs at low cost: GISS 
(Goddard Institute for Space Studies) modelE-R. The atmo- 
sphere model is the documented modelE (Schmidt et al., 
2006). The ocean is based on the Russell et al. (1995) model 
that conserves water and salt mass; has a free surface with 
divergent flow; uses a linear upstream scheme for advection; 
allows flow through 12 sub-resolution straits; and has back- 
ground diffusivity of 0.3 cm 2 s -1 , resolution of 4° x 5° and 
13 layers that increase in thickness with depth. 

However, the ocean model includes simple but significant 
changes, compared with the version documented in simula- 
tions by Miller et al. (2014). First, an error in the calculation 
of neutral surfaces in the Gent-McWilliams (GM; Gent and 
McWilliams, 1990) mesoscale eddy parameterization was 
corrected; the resulting increased slope of neutral surfaces 
provides proper leverage to the restratification process and 
correctly orients eddy stirring along those surfaces. 

Second, the calculation of eddy diffusivity K meso for GM 
following Visbeck et al. (1997) was simplified to use a length 
scale independent of the density structure (J. Marshall, per- 
sonal communication, 2012): 

K meso = C/[7Eady x /(latitude)], (2) 

where C = (27.9km) 2 , Eady growth rate l/7Eady = {|S x 
/V | } , S is the neutral surface slope, N the Brunt-Vaisala 
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frequency, { } signifies averaging over the upper D me- 
ters of ocean depth, D — min(max(depth, 400 m), 1000 m), 
and /(latitude) = max(0.1, sin( | latitude I)) 1 to qualitatively 
mimic the larger values of the Rossby radius of deformation 
at low latitudes. These choices for K meso , whose simplicity is 
congruent with the use of a depth-independent eddy diffusiv- 
ity and the use of 1 /T e a dy as a metric of eddy energy, result in 
the zonal average diffusivity shown in Fig. 1 . Third, the so- 
called nonlocal terms in the KPP mixing parameterization 
(Large et al., 1994) were activated. All of these modifica- 
tions tend to increase the ocean stratification, and in particu- 
lar the Southern Ocean state is fundamentally improved. For 
example, we show in Sect. 3.8.5 that our current model pro- 
duces Antarctic Bottom Water on the Antarctic coastline, as 
observed, rather than in the middle of the Southern Ocean as 
occurs in many models, including the GISS-ER model docu- 
mented in CMIP5. However, although overall realism of the 
ocean circulation is much improved, significant model defi- 
ciencies remain, as we will describe. 

The simulated Atlantic meridional overturning circulation 
(AMOC) has maximum flux that varies within the range 
~ 14-18 Sv in the model control run (Figs. 2 and 3). AMOC 
strength in recent observations is 17.5 ± 1.6 Sv (Baringer et 
al., 2013; Srokosz et al., 2012), based on 8 years (2004- 
2011) of data for an in situ mooring array (Rayner et al., 
2011; Johns etal., 2011). 

Ocean model control run initial conditions are climatology 
for temperature and salinity (Levitus and Boyer, 1994; Lev- 
itus et al., 1994); atmospheric composition is that of 1880 
(Hansen et al., 2011). Overall model drift from control run 
initial conditions is moderate (see Fig. SI for planetary en- 
ergy imbalance and global temperature), but there is drift in 
the North Atlantic circulation. The AMOC circulation cell 
initially is confined to the upper 3 km at all latitudes (1st cen- 
tury in Figs. 2 and 3), but by the 5th century the cell reaches 
deeper at high latitudes. 

Atmospheric and surface climate in the present model is 
similar to the documented modelE-R, but because of changes 
to the ocean model we provide several diagnostics in the Sup- 
plement. A notable flaw in the simulated surface climate is 
the unrealistic double precipitation maximum in the tropi- 
cal Pacific (Fig. S2). This double Intertropical Convergence 
Zone (ITCZ) occurs in many models and may be related to 
cloud and radiation biases over the Southern Ocean (Hwang 
and Frierson, 2013) or deficient low level clouds in the trop- 
ical Pacific (de Szoeke and Xie, 2008). Another flaw is un- 
realistic hemispheric sea ice, with too much sea ice in the 


'Where ocean depth exceeds 1000m. these conditions yield 
D = 1000 m, thus excluding any first-order abyssal bathymetric im- 
print on upper ocean eddy energy, consistent with theory and ob- 
servations. The other objective of the stated condition is to limit 
release of potential energy in the few ocean gridboxes with ocean 
depth less than 400 m, because shallow depths limit the ability of 
baroclinic eddies to release potential energy via vertical motion. 
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Figure 1. Zonal-mean mesoscale diffusivity (m 2 s 1 ) versus time in control run. 


Atlantic meridional overturning circulation (Sv) 
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Figure 2. AMOC (Sv) in the 1st, 5th, 10th, 15th and 20th centuries of the control run. 


Northern Hemisphere and too little in the Southern Hemi- 
sphere (Figs. S3 and S4). Excessive Northern Hemisphere 
sea ice might be caused by deficient poleward heat transport 
in the Atlantic Ocean (Fig. S5). However, the AMOC has re- 
alistic strength and Atlantic meridional heat transport is only 
slightly below observations at high latitudes (Fig. S5). Thus 
we suspect that the problem may lie in sea ice parameteriza- 
tions or deficient dynamical transport of ice out of the Arctic. 
The deficient Southern Hemisphere sea ice, at least in part, 
is likely related to excessive poleward (southward) transport 
of heat by the simulated global ocean (Fig. S5), which is re- 
lated to deficient northward transport of heat in the modeled 
Atlantic Ocean (Fig. S5). 

A key characteristic of the model and the real world is 
the response time: how fast does the surface temperature ad- 
just to a climate forcing? ModelE-R response is about 40 % 
in 5 years (Fig. 4) and 60% in 100 years, with the remain- 
der requiring many centuries. Hansen et al. (2011) concluded 
that most ocean models, including modelE-R, mix a surface 
temperature perturbation downward too efficiently and thus 
have a slower surface response than the real world. The basis 
for this conclusion was empirical analysis using climate re- 
sponse functions, with 50, 75 and 90 % response at year 100 
for climate simulations (Hansen et al., 2011). Earth’s mea- 
sured energy imbalance in recent years and global temper- 
ature change in the past century revealed that the response 
function with 75 % response in 100 years provided a much 
better fit with observations than the other choices. Durack 
et al. (2012) compared observations of how rapidly surface 


salinity changes are mixed into the deeper ocean with the 
large number of global models in CMIP3, reaching a similar 
conclusion, that the models mix too rapidly. 

Our present ocean model has a faster response on 10-75- 
year timescales than the old model (Fig. 4), but the change 
is small. Although the response time in our model is similar 
to that in many other ocean models (Hansen et al., 201 1), we 
believe that it is likely slower than the real-world response 
on timescales of a few decades and longer. A too slow sur- 
face response could result from excessive small-scale mix- 
ing. We will argue, after the studies below, that excessive 
mixing likely has other consequences, e.g., causing the ef- 
fect of freshwater stratification on slowing Antarctic Bottom 
Water (AABW) formation and growth of Antarctic sea ice 
cover to occur 1-2 decades later than in the real world. Sim- 
ilarly, excessive mixing probably makes the AMOC in the 
model less sensitive to freshwater forcing than the real-world 
AMOC. 

3.3 Experiment definition: exponentially increasing 
freshwater 

Freshwater injection is 360 Gt year - 1 ( 1 mm sea level) in 
2003-2015, then grows with 5-, 10- or 20-year doubling time 
(Fig. 5) and terminates when global sea level reaches 1 or 
5 m. Doubling times of 10, 20 and 40 years, reaching meter- 
scale sea level rise in 50, 100, and 200 years may be a more 
realistic range of timescales, but 40 years yields little effect 
this century, the time of most interest, so we learn more with 
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Atlantic meridional overturning circulation (Sv): Control run (2000 years) 



Figure 3. Annual mean AMOC (Sv) at 28° N in the model control run. 



Year 


Figure 4. Climate response function, R(t), i.e., the fraction (%) of 
equilibrium surface temperature response for GISS modelE-R based 
on a 2000-year control run (Hansen et al., 2007a). Forcing was in- 
stant CO 2 doubling with fixed ice sheets, vegetation distribution, 
and other long-lived GHGs. 


less computing time using the 5-, 10- and 20-year doubling 
times. Observed ice sheet mass loss doubling rates, although 
records are short, are ~ 10 years (Sect. 5.1). Our sharp cut- 
off of melt aids separation of immediate forcing effects and 
feedbacks. 

We argue that such a rapid increase in meltwater is plausi- 
ble if GHGs keep growing rapidly. Greenland and Antarctica 
have outlet glaciers in canyons with bedrock below sea level 
well back into the ice sheet (Fretwell et al., 2013; Morlighem 
et al., 2014; Pollard et al., 2015). Feedbacks, including ice 
sheet darkening due to surface melt (Hansen et al., 2007b; 
Robinson et al., 2012; Tedesco et al., 2013; Box et al., 2012) 
and lowering and thus warming of the near-coastal ice sheet 
surface, make increasing ice melt likely. Paleoclimate data 
reveal sea level rise of several meters in a century (Fair- 
banks, 1989; Deschamps et al., 2012). Those cases involved 
ice sheets at lower latitudes, but 21st century climate forcing 
is larger and increasing much more rapidly. 


Radiative forcings (Fig. S16a, b) are from Hansen et 
al. (2007c) through 2003 and IPCC scenario A1B for later 
GHGs. A1B is an intermediate IPCC scenario over the cen- 
tury, but on the high side early this century (Fig. 2 of Hansen 
et al., 2007c). We add freshwater to the North Atlantic (ocean 
area within 52-72° N and 15°E-65°N) or Southern Ocean 
(ocean south of 60° S), or equally divided between the 
two oceans. Ice sheet discharge (icebergs plus meltwater) is 
mixed as freshwater with mean temperature —15 °C into the 
top three ocean layers (Fig. S6). 

3.4 Simulated surface temperature and energy balance 

We present surface temperature and planetary energy balance 
first, thus providing a global overview. Then we examine 
changes in ocean circulation and compare results with prior 
studies. 

Temperature change in 2065, 2080 and 2096 for 10-year 
doubling time (Fig. 6) should be thought of as results when 
sea level rise reaches 0.6, 1 .7 and 5 m, because the dates de- 
pend on initial freshwater flux. Actual current freshwater flux 
may be about a factor of 4 higher than assumed in these 
initial runs, as we will discuss, and thus effects may occur 
~ 20 years earlier. A sea level rise of 5 m in a century is about 
the most extreme in the paleo-record (Fairbanks, 1989; De- 
schamps et al., 2012), but the assumed 21st century climate 
forcing is also more rapidly growing than any known natural 
forcing. 

Meltwater injected into the North Atlantic has larger ini- 
tial impact, but Southern Hemisphere ice melt has a greater 
global effect for larger melt as the effectiveness of more melt- 
water in the North Atlantic begins to decline. The global ef- 
fect is large long before sea level rise of 5 m is reached. Melt- 
water reduces global warming about half by the time sea level 
rise reaches 1 .7 m. Cooling due to ice melt more than elimi- 
nates A1B warming in large areas of the globe. 

The large cooling effect of ice melt does not decrease 
much as the ice melting rate varies between doubling times 
of 5, 10 or 20 years (Fig. 7a). In other words, the cumulative 
ice sheet melt, rather than the rate of ice melt, largely deter- 
mines the climate impact for the range of melt rates covered 
by 5-, 10- and 20-year doubling times. Thus if ice sheet loss 
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(a) Freshwater melt flux (b) Resulting sea level rise 



Figure 5. (a) Total freshwater flux added in the North Atlantic and Southern oceans and (b) resulting sea level rise. Solid lines for 1 m sea 
level rise, dotted for 5 m. One sverdrup (Sv) is 10 6 m 3 s _1 , which is ~ 3x 10 4 Gtyear -1 . 
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Figure 6. Surface air temperature (°C) relative to 1880-1920 in (a) 2065, (b) 2080, and (c) 2096. Top row is IPCC scenario A1B. Ice melt 
with 10-year doubling is added in other scenarios. 


occurs even to an extent of 1.7 m sea level rise (Fig. 7b), a 
large impact on climate and climate change is predicted. 

Greater global cooling occurs for freshwater injected into 
the Southern Ocean, but the cooling lasts much longer for 
North Atlantic injection (Fig. 7a). That persistent cooling, 
mainly at Northern Hemisphere middle and high latitudes 


(Fig. S7), is a consequence of the sensitivity, hysteresis ef- 
fects, and long recovery time of the AMOC (Stocker and 
Wright, 1991; Rahmstorf, 1995, and earlier studies refer- 
enced therein). AMOC changes are described below. 

When freshwater injection in the Southern Ocean is halted, 
global temperature jumps back within two decades to the 
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(a) Global surface air temperature 



Figure 7. (a) Surface air temperature (°C) relative to 1880-1920 for several scenarios, (b) Global energy imbalance (W m 2 ) for the same 
scenarios. 


value it would have had without any freshwater addition 
(Fig. 7a). Quick recovery is consistent with the Southern 
Ocean-centric picture of the global overturning circulation 
(Fig. 4; Talley, 2013), as the Southern Ocean meridional 
overturning circulation (SMOC), driven by AABW forma- 
tion, responds to change in the vertical stability of the ocean 
column near Antarctica (Sect. 3.7) and the ocean mixed layer 
and sea ice have limited thermal inertia. 

Cooling from ice melt is largely regional, temporary, and 
does not alleviate concerns about global warming. Southern 
Hemisphere cooling is mainly in uninhabited regions. North- 
ern Hemisphere cooling increases temperature gradients that 
will drive stronger storms (Sect. 3.9). 

Global cooling due to ice melt causes a large increase in 
Earth’s energy imbalance (Fig. 7b), adding about +2 W m 2 , 
which is larger than the imbalance caused by increasing 
GHGs. Thus, although the cold freshwater from ice sheet 
disintegration provides a negative feedback on regional and 
global surface temperature, it increases the planet’s energy 
imbalance, thus providing more energy for ice melt (Hansen, 
2005). This added energy is pumped into the ocean. 

Increased downward energy flux at the top of the atmo- 
sphere is not located in the regions cooled by ice melt. How- 
ever, those regions suffer a large reduction of net incoming 
energy (Fig. 8a). The regional energy reduction is a conse- 
quence of increased cloud cover (Fig. 8b) in response to the 
colder ocean surface. However, the colder ocean surface re- 
duces upward radiative, sensible and latent heat fluxes, thus 
causing a large (~ 50 W m -2 ) increase in energy into the 


North Atlantic and a substantial but smaller flux into the 
Southern Ocean (Fig. 8c). 

Below we conclude that the principal mechanism by which 
this ocean heat increases ice melt is via its effect on ice 
shelves. Discussion requires examination of how the fresh- 
water injections alter the ocean circulation and internal ocean 
temperature. 

3.5 Simulated Atlantic meridional overturning 
circulation (AMOC) 

Broecker’s articulation of likely effects of freshwater out- 
bursts in the North Atlantic on ocean circulation and global 
climate (Broecker, 1990; Broeckeret al., 1990) spurred quan- 
titative studies with idealized ocean models (Stocker and 
Wright, 1991) and global atmosphere-ocean models (Man- 
abe and Stouffer, 1995; Rahmstorf 1995, 1996). Scores of 
modeling studies have since been carried out, many reviewed 
by Barreiro et al. (2008), and observing systems are being 
developed to monitor modern changes in the AMOC (Carton 
and Hakkinen, 201 1). 

Our climate simulations in this section are five-member 
ensembles of runs initiated at 25-year intervals at years 901- 
1001 of the control run. We chose this part of the control 
run because the planet is then in energy balance (Fig. SI), 
although by that time model drift had altered the slow deep- 
ocean circulation. Some model drift away from initial clima- 
tological conditions is inevitable, as all models are imperfect, 
and we carry out the experiments with cognizance of model 
limitations. However, there is strong incentive to seek basic 
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Change in 2080 relative to 1880-1920 

(a) Planetary energy balance (W/m 2 ) (b) Cloud cover (%) (c) Net energy into ground (W/m 2 ) 
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Figure 8. Change in 2080 (mean of 2078-2082), relative to 1880-1920, of annual mean (a) planetary energy balance (W m 2 ), (b) cloud 
cover (%), and (c) net energy into ground (W m -2 ) for the same scenarios as Fig. 6. 


improvements in representation of physical processes to re- 
duce drift in future versions of the model. 

GHGs alone (scenario A1B) slow AMOC by the early 
21st century (Fig. 9), but variability among individual mns 
(Fig. S8) would make definitive detection difficult at present. 
Freshwater injected into the North Atlantic or in both hemi- 
spheres shuts down the AMOC (Fig. 9, right side). GHG 
amounts are fixed after 2100 and ice melt is zero, but after 
two centuries of stable climate forcing the AMOC has not 
recovered to its earlier state. This slow recovery was found 
in the earliest simulations by Manabe and Stouffer (1994) 
and Rahmstorf (1995, 1996). 

Freshwater injection already has a large impact when ice 
melt is a fraction of 1 m of sea level. By the time sea level rise 
reaches 59 cm (2065 in the present scenarios), when fresh- 
water flux is 0.48 Sv, the impact on AMOC is already large, 
consistent with the substantial surface cooling in the North 
Atlantic (Fig. 6). 


3.6 Comparison with prior simulations 

AMOC sensitivity to GHG forcing has been examined ex- 
tensively based on CMIP studies. Schmittner et al. (2005) 
found that AMOC weakened 25 ± 25 % by the end of the 
21st century in 28 simulations of 9 different models forced by 
the A1B emission scenario. Gregory et al. (2005) found 10- 
50 % AMOC weakening in 1 1 models for COt quadrupling 
(1 % year -1 increase for 140 years), with largest decreases in 
models with strong AMOCs. Weaver et al. (2007) found a 
15-31 % AMOC weakening for CO 2 quadrupling in a single 
model for 17 climate states differing in initial GHG amount. 
AMOC in our model weakens 30 % in the century between 
1990-2000 and 2090-2100, the period used by Schmittner 
et al. (2005), for A1B forcing (Fig. S8). Thus our model is 
more sensitive than the average but within the range of other 
models, a conclusion that continues to be valid in comparison 
with 10 CMIP5 models (Cheng et al., 2013). 

AMOC sensitivity to freshwater forcing has not been com- 
pared as systematically among models. Several studies find 
little impact of Greenland melt on AMOC (Huybrechts et 
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Atlantic meridional overturning circulation (Sv) at 28° N 



Ice melt in Southern Ocean 


Ice melt in both hemispheres 


Figure 9. Ensemble-mean AMOC (Sv) at 28° N versus time for the same four scenarios as in Fig. 6, with ice melt reaching 5 m at the end of 
the 21st century in the three experiments with ice melt. 



al., 2002; Jungclaus et al., 2006; Vizcaino et al., 2008) while 
others find substantial North Atlantic cooling (Fichefet et al., 
2003; Swingedouw et al., 2007; Hu et al., 2009, 201 1). Stud- 
ies with little impact calculated or assumed small ice sheet 
melt rates, e.g., Greenland contributed only 4 cm of sea level 
rise in the 21st century in the ice sheet model of Huybrechts 
et al. (2002). Fichefet et al. (2003), using nearly the same 
atmosphere-ocean model as Huybrechts et al. (2002) but a 
more responsive ice sheet model, found AMOC weakening 
from 20 to 13 Sv late in the 21st century, but separate contri- 
butions of ice melt and GHGs to AMOC slowdown were not 
defined. 

Hu et al. (2009, 2011) use the A1B scenario and freshwa- 
ter from Greenland starting at 1 mm sea level per year in- 
creasing 7% year -1 , similar to our 10-year doubling case. 
Hu et al. keep the melt rate constant after it reaches 0.3 Sv 
(in 2050), yielding 1.65 m sea level rise in 2100 and 4.2m in 
2200. Global warming found by Hu et al. for scenario A1B 
resembles our result but is 20-30 % smaller (compare Fig. 2b 
of Hu et al., 2009 to our Fig. 6), and cooling they obtain from 
the freshwater flux is moderately less than that in our model. 
AMOC is slowed about one-third by the latter 21st century 
in the Hu et al. (2011) 7 % year -1 experiment, comparable to 
our result. 

General consistency holds for other quantities, such as 
changes of precipitation. Our model yields southward shift- 
ing of the Intertropical Convergence Zone (ITCZ) and in- 
tensification of the subtropical dry region with increasing 
GHGs (Fig. S9), as has been reported in modeling studies 
of Swingedouw et al. (2007, 2009). These effects are inten- 
sified by ice melt and cooling in the North Atlantic region 
(Fig. S9). 

A recent five-model study (Swingedouw et al., 2014) finds 
a small effect on AMOC for 0. 1 Sv Greenland freshwater flux 
added in 2050 to simulations with a strong GHG forcing. Our 


larger response is likely due, at least in part, to our freshwater 
flux reaching several tenths of a sverdrup. 

3.7 Pure freshwater experiments 

We assumed, in discussing the relevance of these experi- 
ments to Eemian climate, that effects of freshwater injec- 
tion dominate over changing GHG amount, as seems likely 
because of the large freshwater effect on sea surface tem- 
peratures (SSTs) and sea level pressure. However, Eemian 
CO 2 was actually almost constant at ~ 275 ppm (Luthi et al., 
2008). Thus, to isolate effects better, we now carry out sim- 
ulations with fixed GHG amount, which helps clarify impor- 
tant feedback processes. 

Our pure freshwater experiments are five-member ensem- 
bles starting at year's 1001, 1101, 1201, 1301, and 1401 of the 
control run. Each experiment ran 300 years. Freshwater flux 
in the initial decade averaged 180 km 3 year -1 (0.5 mm sea 
level) in the hemisphere with ice melt and increased with a 
10-year doubling time. Freshwater input is terminated when 
it reaches 0.5 m sea level rise per hemisphere for three five- 
member ensembles: two ensembles with injection in the in- 
dividual hemispheres and one ensemble with input in both 
hemispheres (lm total sea level rise). Three additional en- 
sembles were obtained by continuing freshwater injection 
until hemispheric sea level contributions reached 2.5 m. Here 
we provide a few model diagnostics central to discussions 
that follow. Additional results are provided in Figs. S 10-S 12. 

The AMOC shuts down for Northern Hemisphere fresh- 
water input yielding 2.5 m sea level rise (Fig. 10). By year 
300, more than 200 years after cessation of all freshwater in- 
put, AMOC is still far from full recovery for this large fresh- 
water input. On the other hand, freshwater input of 0.5 m 
does not cause full shutdown, and AMOC recovery occurs 
in less than a century. 
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Atlantic meridional overturning circulation (Sv) at 28° N 

Sea level rise to Sea level rise to 





Figure 10. Ensemble-mean AMOC (Sv) at 28° N versus time for six pure freshwater forcing experiments. 


Global temperature change (Fig. 11) reflects the funda- 
mentally different impact of freshwater forcings of 0.5 and 
2.5 m. The response also differs greatly depending on the 
hemisphere of the freshwater input. The case with freshwater 
forcing in both hemispheres is shown only in the Supplement 
because, to a good approximation, the response is simply the 
sum of the responses to the individual hemispheric forcings 
(see Figs. S10-S12). The sum of responses to hemispheric 
forcings moderately exceeds the response to global forcing. 

Global cooling continues for centuries for the case with 
freshwater forcing sufficient to shut down the AMOC 
(Fig. 11). If the forcing is only 0.5 m of sea level, the tem- 
perature recovers in a few decades. However, the freshwater 
forcing required to reach the tipping point of AMOC shut- 
down may be less in the real world than in our model, as 
discussed below. Global cooling due to freshwater input on 
the Southern Ocean disappears in a few years after freshwa- 
ter input ceases (Fig. 1 1), for both the smaller (0.5 m of sea 
level) and larger (2.5 m) freshwater forcings. 

Injection of a large amount of surface freshwater in ei- 
ther hemisphere has a notable impact on heat uptake by the 
ocean and the internal ocean heat distribution (Fig. 12). De- 
spite continuous injection of a large amount of very cold 
(— 15 °C) water in these pure freshwater experiments, sub- 
stantial portions of the ocean interior become warmer. Trop- 
ical and Southern Hemisphere warming is the well-known 
effect of reduced heat transport to northern latitudes in re- 
sponse to the AMOC shutdown (Rahmstorf, 1996; Barreiro 
et al., 2008). 


However, deep warming in the Southern Ocean may have 
greater consequences. Warming is maximum at grounding 
line depths (~ 1-2 km) of Antarctic ice shelves (Rignot and 
Jacobs, 2002). Ice shelves near their grounding lines (Fig. 13 
of Jenkins and Doake, 1991) are sensitive to temperature of 
the proximate ocean, with ice shelf melting increasing 1 m 
per year for each 0. 1 °C temperature increase (Rignot and Ja- 
cobs, 2002). The foot of an ice shelf provides most of the re- 
straining force that ice shelves exert on landward ice (Fig. 14 
of Jenkins and Doake, 1991), making ice near the grounding 
line the buttress of the buttress. Pritchard et al. (2012) deduce 
from satellite altimetry that ice shelf melt has primary control 
of Antarctic ice sheet mass loss. 

Thus we examine our simulations in more detail (Fig. 13). 
The pure freshwater experiments add 5 mm sea level in the 
first decade (requiring an initial 0.346 mm year -1 for 10-year 
doubling), 10 mm in the second decade, and so on (Fig. 13a). 
Cumulative freshwater injection reaches 0.5 m in year 68 and 
2.5 m in year 90. 

Antarctic Bottom Water (AABW) formation is reduced 
~ 20 % by year 68 and ~ 50 % by year 90 (Fig. 13b). When 
freshwater injection ceases, AABW formation rapidly re- 
gains full strength, in contrast to the long delay in reestab- 
lishing North Atlantic Deep Water (NADW) formation af- 
ter AMOC shutdown. The Southern Ocean mixed-layer re- 
sponse time dictates the recovery time for AABW forma- 
tion. Thus rapid recovery also applies to ocean temperature 
at depths of ice shelf grounding lines (Fig. 13c). The rapid 
response of the Southern Ocean meridional overturning cir- 
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Global surface air temperature 


(a) Freshwater in North Atlantic (b) Freshwater in Southern Ocean 



Figure 11. Ensemble-mean global surface air temperature (°C) for experiments (years on x axis) with freshwater forcing in either the North 
Atlantic Ocean (left) or the Southern Ocean (right). 
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Figure 12. Change of ocean temperature (°C) relative to control run due to freshwater input that reaches 2.5 m of global sea level in a 
hemisphere (thus 5 m sea level rise in the bottom row). 


culation (SMOC) implies that the rate of freshwater addition 
to the mixed layer is the driving factor. 

Freshwater flux has little effect on simulated Northern 
Hemisphere sea ice until the 7th decade of freshwater growth 
(Fig. 13d), but Southern Hemisphere sea ice is more sensi- 
tive, with substantial response in the 5th decade and large re- 
sponse in the 6th decade. Below we show that “5th decade” 
freshwater flux (2880 Cit year’ ' ) is already relevant to the 
Southern Ocean today. 


3.8 Simulations to 2100 with modified (more realistic) 
forcings 

Recent data show that current ice melt is larger than assumed 
in our 1850-2300 simulations. Thus we make one more sim- 
ulation and include minor improvements in the radiative forc- 
ing. 

3.8.1 Advanced (earlier) freshwater injection 

Atmosphere-ocean climate models, including ours, com- 
monly include a fixed freshwater flux from the Greenland 
and Antarctic ice sheets to the ocean. This flux is chosen to 
balance snow accumulation in the model’s control run, with 
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(a) Rate of freshwater input (b) SMOC (mean ocean overturning at 72°S) 




Figure 13. (a) Freshwater input (Tt year 1 ) to Southern Ocean (1 Tt = 1000 km 3 ), (b, c, d) Simulated overturning strength (Sv) of AABW 
cell at 72° S, temperature (°C) at depth 1 . 13 km at 74° S, and sea ice cover (%). 


the rationale that approximate balance is expected between 
net accumulation and mass loss including icebergs and ice 
shelf melting. Global warming creates a mass imbalance that 
we want to investigate. Ice sheet models can calculate the im- 
balance, but it is unclear how reliably ice sheet models sim- 
ulate ice sheet disintegration. We forgo ice sheet modeling, 
instead adding a growing freshwater amount to polar oceans 
with alternative growth rates and initial freshwater amount 
estimated from available data. 

Change of freshwater flux into the ocean in a warming 
world with shrinking ice sheets consists of two terms: term 
1 being net ice melt and term 2 being change in P — E (pre- 
cipitation minus evaporation) over the relevant ocean. Term 
1 includes land based ice mass loss, which can be detected 
by satellite gravity measurements, loss of ice shelves, and 
net sea ice mass change. Term 2 is calculated in a climate 
model forced by changing atmospheric composition, but it is 
not included in our pure freshwater experiments that have no 
global warming. 

IPCC (Vaughan et al., 2013) estimated land ice loss in 
Antarctica that increased from 30Gtyear -1 in 1992-2001 
to 1 47 Gt year - 1 in 2002-2011 and in Greenland from 34 
to 215 Gt year -1 , with uncertainties discussed by Vaughan 
et al. (2013). Gravity satellite data suggest Greenland ice 
sheet mass loss ~ 300-400 Gt year - 1 in the past few years 
(Barletta et al., 2013). A newer analysis of gravity data for 
2003-2013 (Velicogna et al., 2014), discussed in more detail 
in Sect. 5.1, finds a Greenland mass loss 280 ± 58 Gt year -1 
and Antarctic mass loss 67 ± 44 Gt year - 1 . 

One estimate of net ice loss from Antarctica, including 
ice shelves, is obtained by surveying and adding the mass 
flux from all ice shelves and comparing this freshwater mass 
loss with the freshwater mass gain from the continental sur- 
face mass budget. Rignot et al. (2013) and Depoorter et 
al. (2013) independently assessed the freshwater mass fluxes 
from Antarctic ice shelves. Their respective estimates for 
the basal melt are 1500±237 and 1454± 174 Gt year -1 . 


Their respective estimates for calving are 1265 ±139 and 
1321 ± 144 Gt year -1 . 

This estimated freshwater loss via the ice shelves 
(~ 2800 Gt year -1 ) is larger than freshwater gain by Antarc- 
tica. Vaughan et al. (1999) estimated net surface mass bal- 
ance of the continent as ±1811 and ±2288 Gt year -1 in- 
cluding precipitation on ice shelves. Vaughan et al. (2013) 
estimates the net Antarctic surface mass balance as 
±1983 ± 122 Gtyear -1 excluding ice shelves. Thus compar- 
ison of continental freshwater input with ice shelf output sug- 
gests a net export of freshwater to the Southern Ocean of sev- 
eral hundred Gtyear -1 in recent years. However, substantial 
uncertainty exists in the difference between these two large 
numbers. 

An independent evaluation has recently been achieved by 
Rye et al. (2014) using satellite measured changes of sea 
level around Antarctica in the period 1992-2011. Sea level 
along the Antarctic coast rose 2 mm year -1 faster than the 
regional mean sea level rise in the Southern Ocean south of 
50° S, an effect that they conclude is almost entirely a steric 
adjustment caused by accelerating freshwater discharge from 
Antarctica. They conclude that an excess freshwater input of 
430 ±230 Gtyear -1 , above the rate needed to maintain a 
steady ocean salinity, is required. Rye et al. (2014) note that 
these values constitute a lower bound for the actual excess 
discharge above a “steady salinity” rate, because numerous 
in situ data, discussed below, indicate that freshening began 
earlier than 1992. 

Term 2, change in P — E over the Southern Ocean rela- 
tive to its preindustrial amount, is large in our climate sim- 
ulations. In our ensemble of runs (using observed GHGs 
for 1850-2003 and scenario A1B thereafter) the increase in 
P — E in the decade 2011-2020, relative to the control run, 
was in the range 3500 to 4000 Gtyear -1 , as mean precipita- 
tion over the Southern Ocean increased ~ 35 mm year -1 and 
evaporation decreased ~ 3 mm year - 1 . 

Increasing ice melt and increasing P — E are climate feed- 
backs, their growth in recent decades driven by global warm- 
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ing. Our pure freshwater simulations indicate that their sum, 
at least 4000 Gt year - 1 , is sufficient to affect ocean circula- 
tion, sea ice cover, and surface temperature, which can spur 
other climate feedbacks. We investigate these feedbacks via 
climate simulations using improved estimates of freshwater 
flux from ice melt. P — E is computed by the model. 

We take freshwater injection to be 720 Gt year -1 from 
Antarctica and 360 Gt year -1 in the North Atlantic in 2011, 
with injection rates at earlier and later times defined by 
assumption of a 10-year doubling time. Resulting mean 
freshwater injection around Antarctica in 1992-2011 is 
~400Gtyear -1 , similar to the estimate of Rye et al. (2014). 
A recent estimate of 310 ± 74km 3 volume loss of floating 
Antarctic ice shelves in 2003-2012 (Paolo et al., 2015) is 
not inconsistent, as the radar altimeter data employed for ice 
shelves do not include contributions from the ice sheet or 
fast ice tongues at the ice shelf grounding line. Greenland ice 
sheet mass loss provides most of the assumed 360 Gt year -1 
freshwater, and this would be supplemented by shrinking ice 
shelves (Rignot and Steffen, 2008) and small ice caps in the 
North Atlantic and west of Greenland (Ohmura, 2009) that 
are losing mass (Abdalati et al., 2004; Bahr et al., 2009). 

We add freshwater around Antarctica at coastal grid boxes 
(Fig. S13) guided by the data of Rignot et al. (2013) and De- 
poorter et al. (2013). Injection in the Western Hemisphere, 
especially from the Weddell Sea to the Ross Sea, is more 
than twice that in the other hemisphere (Fig. 14). Specified 
freshwater flux around Greenland is similar on the east and 
west coasts, and small along the north coast (Fig. SI 3). 

3.8.2 Modified radiative forcings 

Actual GHG forcing is less than scenario A1B, because CH 4 
and minor gas growth declined after IPCC scenarios were 
defined (Fig. 5; Hansen et al., 2013c, update at http://www. 
columbia.edu/~mhsll9/GHGs/). As a simple improvement 
we decreased the A1B CH 4 scenario during 2003-2013 so 
that subsequent CH 4 is reduced lOOppb, decreasing radiative 
forcing ~ 0.05 W m -2 . 

Stratospheric aerosol forcing to 2014 uses the data set of 
Sato et al. (1993) as updated at http://www.columbia.edu/ 
-rnhsl 19/StratAer/. Future years have constant aerosol op- 
tical depth 0.0052 yielding effective forcing — 0.12Wm -2 , 
implemented by using fixed 1997 aerosol data. Tropospheric 
aerosol growth is assumed to slow smoothly, leveling out at 
— 2 W m -2 in 2100. Future solar forcing is assumed to have 
an 1 1-year cycle with amplitude 0.25 W m -2 . Net forcing ex- 
ceeds 5 W m -2 by the end of the 21st century, about 3 times 
the current forcing (Fig. S16). 

3.8.3 Climate simulations with modified forcings 

Global temperature has a maximum at +1.2 °C in the 2040s 
for the modified forcings (Fig. 15). Ice melt cooling is ad- 
vanced as global ice melt reaches 1 m of sea level in 2060, 
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Figure 14. Freshwater flux (Gt year - 1 ) from Antarctic ice shelves 
based on data of Rignot et al. (2013), integrated here into intervals 
of 15° of longitude. Depoorter et al. (2013) data yield a similar 
distribution. 


1/3 from Greenland and 2/3 from Antarctica. Global tem- 
perature rise resumes in the 2060s after cessation of freshwa- 
ter injection. 

Global temperature becomes an unreliable diagnostic of 
planetary condition as the ice melt rate increases. Global en- 
ergy imbalance (Fig. 15b) is a more meaningful measure of 
planetary status as well as an estimate of the climate forcing 
change required to stabilize climate. Our calculated present 
energy imbalance of ~ 0.8 W m -2 (Fig. 15b) is larger than 
the observed 0.58 ± 0.15 W m -2 during 2005-2010 (Hansen 
et al., 2011). The discrepancy is likely accounted for by ex- 
cessive ocean heat uptake at low latitudes in our model, a 
problem related to the model’s slow surface response time 
(Fig. 4) that may be caused by excessive small-scale ocean 
mixing. 

Large scale regional cooling occurs in the North Atlantic 
and Southern oceans by mid-century (Fig. 16) for 10-year 
doubling of freshwater injection. A 20-year doubling places 
similar cooling near the end of this century, 40 years ear- 
lier than in our prior simulations (Fig. 7), as the factor of 4 
increase in current freshwater from Antarctica is a 40-year 
advance. 

Cumulative North Atlantic freshwater forcing in sverdrup 
years (Sv years) is 0.2 Sv years in 2014, 2.4 Sv years in 2050, 
and 3.4 Sv years (its maximum) prior to 2060 (Fig. S14). 
The critical issue is whether human-spurred ice sheet mass 
loss can be approximated as an exponential process during 
the next few decades. Such nonlinear behavior depends upon 
amplifying feedbacks, which, indeed, our climate simula- 
tions reveal in the Southern Ocean. 

3.8.4 Southern Ocean feedbacks 

Amplifying feedbacks in the Southern Ocean and atmo- 
sphere contribute to dramatic climate change in our simula- 
tions (Fig. 16). We first summarize the feedbacks to identify 
processes that must be simulated well to draw valid conclu- 
sions. While recognizing the complexity of the global ocean 
circulation (Lozier, 2012; Lumpkin and Speer, 2007; Mar- 
shall and Speer, 2012; Munk and Wunsch, 1998; Orsi et al., 
1999; Sheen et al., 2014; Talley, 2013; Wunsch and Ferrari, 
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(a) Global surface air temperature 



Figure 15. (a) Surface air temperature (°C) change relative to 1880-1920 and (b) global energy imbalance (W m ~) for the modified forcing 
scenario including cases with global ice melt reaching 1 and 3.8 m. 


2055-2060 surface air temperature (°C) relative to 1880-1920 


A1B + modified forcings, ice melt to 1 m 1.19 



Figure 16. Surface air temperature (°C) change relative to 1880— 
1920 in 2055-2060 for modified forcings. 


2004), we use a simple two-dimensional representation to 
discuss the feedbacks. 

Climate change includes slowdown of AABW formation, 
indeed shutdown by mid-century if freshwater injection in- 
creases with a doubling time as short as 10 years (Fig. 17). 
Implications of AABW shutdown are so great that we must 
ask whether the mechanisms are simulated with sufficient re- 
alism in our climate model, which has coarse resolution and 
relevant deficiencies that we have noted. After discussing the 
feedbacks here, we examine how well the processes are in- 
cluded in our model (Sect. 3.8.5). Paleoclimate data (Sect. 4) 


provide much insight about these processes, and modern ob- 
servations (Sect. 5) suggest that these feedbacks are already 
underway. 

Large-scale climate processes affecting ice sheets are 
sketched in Fig. 18. The role of the ocean circulation in the 
global energy and carbon cycles is captured to a useful extent 
by the two-dimensional (zonal-mean) overturning circulation 
featuring deep water (NADW) and bottom water (AABW) 
formation in the polar regions. Marshall and Speer (2012) 
discuss the circulation based in part on tracer data and analy- 
ses by Lumpkin and Speer (2007). Talley (2013) extends the 
discussion with diagrams clarifying the role of the Pacific 
and Indian oceans. 

Wunsch (2002) emphasizes that the ocean circulation is 
driven primarily by atmospheric winds and secondarily by 
tidal stirring. Strong circumpolar westerly winds provide en- 
ergy drawing deep water toward the surface in the South- 
ern Ocean. Ocean circulation also depends on processes 
maintaining the ocean’s vertical density stratification. Winter 
cooling of the North Atlantic surface produces water dense 
enough to sink (Fig. 18), forming North Atlantic Deep Wa- 
ter (NADW). However, because North Atlantic water is rel- 
atively fresh, compared to the average ocean, NADW does 
not sink all the way to the global ocean bottom. Bottom 
water is formed instead in the winter around the Antarctic 
coast, where very salty cold water (AABW) can sink to the 
ocean floor. This ocean circulation (Fig. 18) is altered by nat- 
ural and human-made forcings, including freshwater from ice 
sheets, engendering powerful feedback processes. 
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Global meridional overturning circulation (Sv) at 72° S 



Figure 17. SMOC, ocean overturning strength (Sv) at 72° S, including only the mean (Eulerian) transport. This is the average of a five- 
member model ensemble for the modified forcing including advanced ice melt (720 Gt year -1 from Antarctica in 2011) and 10-year doubling. 



Figure 18. Schematic of stratification and precipitation amplifying feedbacks. Stratification: increased freshwater flux reduces surface water 
density, thus reducing AABW formation, trapping NADW heat, and increasing ice shelf melt. Precipitation: increased freshwater flux cools 
ocean mixed layer, increases sea ice area, causing precipitation to fall before it reaches Antarctica, reducing ice sheet growth and increasing 
ocean surface freshening. Ice in West Antarctica and the Wilkes Basin, East Antarctica, is most vulnerable because of the instability of 
retrograde beds. 


A key Southern Ocean feedback is meltwater stratification 
effect, which reduces ventilation of ocean heat to the atmo- 
sphere and space. Our “pure freshwater” experiments show 
that the low-density lid causes deep-ocean warming, espe- 
cially at depths of ice shelf grounding lines that provide most 
of the restraining force limiting ice sheet discharge (Fig. 14 
of Jenkins and Doake, 1991). West Antarctica and Wilkes 
Basin in East Antarctica have potential to cause rapid sea 
level rise, because much of their ice sits on retrograde beds 
(beds sloping inland), a situation that can lead to unstable 
grounding line retreat and ice sheet disintegration (Mercer, 
1978). 

Another feedback occurs via the effect of surface and 
atmospheric cooling on precipitation and evaporation over 
the Southern Ocean. CMIP5 climate simulations, which do 


not include increasing freshwater injection in the Southern 
Ocean, find snowfall increases on Antarctica in the 21st cen- 
tury, thus providing a negative term to sea level change. 
Frieler et al. (2015) note that 35 climate models are consis- 
tent in showing that warming climate yields increasing snow 
accumulation in accord with paleo-data for warmer climates, 
but the paleo-data refer to slowly changing climate in quasi- 
equilibrium with ocean boundary conditions. In our experi- 
ments with growing freshwater injection, the increasing sea 
ice cover and cooling of the Southern Ocean surface and at- 
mosphere cause the increased precipitation to occur over the 
Southern Ocean, rather than over Antarctica. This feedback 
not only reduces any increase in snowfall over Antarctica but 
also provides a large freshening term to the surface of the 
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Maximum winter mixed layer depth (km) 


Maximum winter mixed layer depth (%) 


February in NH and August in SH 


0.18 February in NF1 and August in SH 


9.50 




Figure 19. Maximum mixed-layer depth (in km, left, and % of ocean depth, right) in February (Northern Hemisphere) and August (Southern 
Hemisphere) using the mixed-layer definition of Heuze et al. (2013). 


Southern Ocean, thus magnifying the direct freshening effect 
from increasing ice sheet melt. 

North Atlantic meltwater stratification effects are also im- 
portant, but different. Meltwater from Greenland can slow 
or shutdown NADW formation, cooling the North Atlantic, 
with global impacts even in the Southern Ocean, as we will 
discuss later. One important difference is that the North At- 
lantic can take centuries to recover from NADW shutdown, 
while the Southern Ocean recovers within 1-2 decades after 
freshwater injection stops (Sect. 3.7). 

3.8.5 Model’s ability to simulate these feedbacks 

Realistic representation of these feedbacks places require- 
ments on both the atmosphere and ocean components of 
our climate model. We discuss first the atmosphere, then the 
ocean. 

There are two main requirements on the atmospheric 
model. First, it must simulate P — E well, because of its im- 
portance for ocean circulation and the amplifying feedback 
in the Southern Ocean. Second, it must simulate winds well, 
because these drive the ocean. 

Simulated P — E (Fig. S 15b) agrees well with meteorolog- 
ical reanalysis (Fig. 3.4b of Rhein et al., 2013). Resulting sea 
surface salinity (SSS) patterns in the model (Fig. S 15a) agree 
well with global ocean surface salinity patterns (Antonov et 
al., 2010, and Fig. 3.4a of Rhein et al., 2013). SSS trends in 
our simulation (Fig. S15c), with the Pacific on average be- 
coming fresher while most of the Atlantic and the subtropics 
in the Southern Hemisphere become saltier, are consistent 
with observed salinity trends (Durack and Wijffels, 2010). 
Recent freshening of the Southern Ocean in our simulation 
is somewhat less than in observed data (Fig. 3.4c, d of Rhein 
et al., 2013), implying that the amplifying feedback may be 
imffe /'estimated in our simulation. A likely reason for that is 
discussed below in conjunction with observed sea ice change. 

Obtaining accurate winds requires the model to simulate 
well atmospheric pressure patterns and their change in re- 


sponse to climate forcings. A test is provided by observed 
changes of the Southern Annular Mode (SAM), with a de- 
crease in surface pressure near Antarctica and a small in- 
crease at midlatitudes (Marshall, 2003) that Thompson et 
al. (2011) relate to stratospheric ozone loss and increasing 
GHGs. Our climate forcing (Fig. S16) includes ozone change 
(Fig. 2 of Hansen et al., 2007a) with stratospheric ozone 
depletion in 1979-1997 and constant ozone thereafter. Our 
model produces a trend toward the high index polarity of 
SAM (Fig. S17) similar to observations, although perhaps 
a slightly smaller change than observed (compare Fig. S17 
with Fig. 3 of Marshall, 2003). SAM continues to increase 
in our model after ozone stabilizes (Fig. S17), suggesting 
that GHGs may provide a larger portion of the SAM re- 
sponse in our model than in the model study of Thompson et 
al. (2011). It would not be surprising if the stratospheric dy- 
namical response to ozone change were weak in our model, 
given the coarse resolution and simplified representation of 
atmospheric drag and dynamical effects in the stratosphere 
(Hansen et al., 2007a), but that is not a major concern for our 
present purposes. 

The ocean model must be able to simulate realistically the 
ocean’s overturning circulation and its response to forcings 
including freshwater additions. Heuze et al. (2013, 2015) 
point out that simulated deep convection in the Southern 
Ocean is unrealistic in most models, with AABW formation 
occurring in the open ocean where it rarely occurs in nature. 
Our present ocean model contains significant improvements 
(see Sect. 3.2) compared to the GISS E2-R model that Heuze 
et al. include in their comparisons. Thus we show (Fig. 19) 
the maximum mixed-layer depth in winter (February in the 
Northern Hemisphere and August in the Southern Hemi- 
sphere) using the same criterion as Heuze et al. to define the 
mixed-layer depth, i.e., the layers with a density difference 
from the ocean surface layer less than 0.03 kg m -3 . 

Southern Ocean mixing in the model reaches a depth of 
~ 500 m in a wide belt near 60° S stretching west from the 
southern tip of South America, with similar depths south 
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of Australia. These open-ocean mixed-layer depths compare 
favorably with observations shown in Fig. 2a of Heuze et 
al. (2015), based on data of de Boyer Montegut et al. (2004). 
There is no open-ocean deep convection in our model. 

Deep convection occurs only along the coast of Antarc- 
tica (Fig. 19). Coastal grid boxes on the continental shelf are 
a realistic location for AABW formation. Orsi et al. (1999) 
suggest that most AABW is formed on shelves around the 
Weddell-Enderby Basin (60 %) and shelves of the Adelie- 
Wilkes Coast and Ross Sea (40 %). Our model produces mix- 
ing down to the shelf in those locations (Fig. 19b), and also 
on the Amery Ice Shelf near the location where Ohshima 
et al. (2013) identified AABW production, which they term 
Cape Darnley Bottom Water. 

With our coarse 4° stair step to the ocean bottom, AABW 
cannot readily slide down the slope to the ocean floor. Thus 
dense shelf water mixes into the open-ocean grid boxes, 
making our modeled Southern Ocean less stratified than the 
real world (cf. temporal drift of Southern Ocean salinity 
in Fig. SI 8), because the denser water must move several 
degrees of latitude horizontally before it can move deeper. 
Nevertheless, our Southern Ocean is sufficiently stratified 
to avoid the unrealistic open-ocean convection that infects 
many models (Heuze et al., 2013, 2015). 

Orsi et al. (1999) estimate the AABW formation rate in 
several ways, obtaining values in the range 8-12 Sv, larger 
than our modeled 5-6 Sv (Fig. 17). However, as in most mod- 
els (Heuze et al., 2015), our SMOC diagnostic (Fig. 17) is 
the mean (Eulerian) circulation, i.e., excluding eddy-induced 
transport. Rerun of a 20-year segment of our control run to 
save eddy-induced changes reveals an increase in SMOC at 
72° S by 1-2 Sv, with negligible change at middle and low 
latitudes, making our simulated transport close to the range 
estimated by Orsi et al. (1999). 

We conclude that the model may simulate Southern Ocean 
feedbacks that magnify the effect of freshwater injected into 
the Southern Ocean: the P — E feedback that wrings global- 
warming-enhanced water vapor from the air before it reaches 
Antarctica and the AABW slowdown that traps deep-ocean 
heat, leaving that heat at levels where it accelerates ice shelf 
melting. Indeed, we will argue that both of these feedbacks 
are probably underestimated in our current model. 

The model seems less capable in Northern Hemisphere 
polar regions. Deep convection today is believed to occur 
mainly in the Greenland-Iceland-Norwegian (GIN) seas and 
at the southern end of Baffin Bay (Fig. 2b of Heuze et al., 
2015). In our model, perhaps because of excessive sea ice 
in those regions, open-ocean deep convection occurs to the 
southeast of the southern tip of Greenland and at less deep 
grid boxes between that location and the United Kingdom 
(Fig. 19). Mixing reaching the ocean floor on the Siberian 
coast in our model (Fig. 19) may be realistic, as coastal 
polynya are observed on the Siberian continental shelf (D. 
Bauch et al., 2012). However, the winter mixed layer on the 
Alaska south coast is unrealistically deep (Fig. 19). These 
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model limitations must be kept in mind in interpreting simu- 
lated Northern Hemisphere climate change. 

3.9 Impact of ice melt on storms 

Our inferences about potential storm changes from continued 
high growth of atmospheric GHGs are fundamentally differ- 
ent than modeling results described in IPCC (2013, 2014), 
where the latter are based on CMIP5 climate model results 
without substantial ice sheet melt. Lehmann et al. (2014) note 
ambiguous results for storm changes from prior model stud- 
ies and describe implications of the CMIP5 ensemble of cou- 
pled climate models. Storm changes are moderate in nature, 
with even a weakening of storms in some locations and sea- 
sons. This is not surprising, because warming is greater at 
high latitudes, reducing meridional temperature gradients. 

Before describing our model results, we note the model 
limitations for study of storms, including its coarse reso- 
lution (4° x 5°), which may contribute to slight misplace- 
ment of the Bermuda high-pressure system for today’s cli- 
mate (Fig. S2). Excessive Northern Hemisphere sea ice may 
cause a bias in location of deepwater formation toward lower 
latitudes. Simulated effects also depend on the location cho- 
sen for freshwater injection; in model results shown here 
(Fig. 20), freshwater was spread uniformly over all longi- 
tudes in the North Atlantic between 65° W and 15° E. It 
would be useful to carry out similar studies with higher- 
resolution models including the most realistic possible dis- 
tribution of meltwater. 

Despite these caveats, we have shown that the model real- 
istically simulates meridional changes of sea level pressure 
in response to climate forcings (Sect. 3.8.5). Specifically, the 
model yields a realistic trend to the positive phase of the 
Southern Annular Mode (SAM) in response to a decrease in 
stratospheric ozone and increase in other GHGs (Fig. S17). 
We also note that the modeled response of atmospheric pres- 
sure to the cooling effect of ice melt is large scale, tending 
to be of a meridional nature that should be handled by our 
model resolution. 

Today’s climate, not Eemian climate, is the base climate 
state upon which we inject polar freshwater. However, the 
simulated climate effects of the freshwater are so large that 
they should also be relevant to freshwater injection in the 
Eemian period. 

3.9.1 Modeling insights into Eemian storms 

Ice melt in the North Atlantic increases simulated sea level 
pressure in that region in all seasons (Fig. 20). In summer 
the Bermuda high-pressure system (Fig. S2) increases in 
strength and moves northward. Circulation around the high 
pressure creates stronger prevailing northeasterly winds at 
latitudes of Bermuda and the Bahamas. A1B climate forc- 
ing alone (Fig. S21, top row) has only a small impact on 
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Sea level pressure and wind: Ice melt in both hemispheres [2078-2082] 

(a) Sea level pressure change (hPa) (b) Wind speed change (m/s) (c) Wind speed (m/s) 


Dec-Jan-Feb 


-0.17 Dec-Jan-Feb 


0.10 Dec-Jan-Feb 


5.35 



Figure 20. Change of seasonal-mean (a) sea level pressure (hPa), (b) wind speed (ms 1 ) in 2078-2082 relative to 1880-1920, and (c) the 
wind speed (ms -1 ) itself, all for the scenario with ice melt in both hemispheres. 


the winds, but cold meltwater in the North Atlantic causes 
a strengthening and poleward shift of the high pressure. 

The high pressure in the model is located further east than 
needed to produce the fastest possible winds at the Bahamas. 
Our coarse-resolution (4° x 5°) model may be partly respon- 
sible for the displacement. However, the location of high 
pressure also depends on meltwater placement, which we 
spread uniformly over all longitudes in the North Atlantic be- 
tween 65° W and 15° E, and on the specific location of ocean 
currents and surface temperature during the Eemian. 

North Atlantic cooling from AMOC shutdown creates 
faster winds in our simulations, with a seasonal-mean incre- 
ment as much as 10-20%. Such a percentage translates into 
an increase in storm power dissipation by a factor ~ 1 .4-2, 
because dissipation is proportional to the cube of wind speed 
(Emanuel, 1987, 2005). Our simulated changes refer to mean 
winds over large grid boxes, not individual storms, for which 
the change in the most extreme cases might be larger. 

Increased North Atlantic high pressure strengthens pre- 
vailing northeasterly winds blowing onto the Bahamas in 


the direction of Eemian wave-formed deposits (Sect. 4.1.2). 
Consistent increase in these winds would contribute to cre- 
ation of long-wavelength, deep-ocean waves that scour the 
ocean floor as they reach the shallow near-shore region. 
However, extreme events may require the combined effect of 
increased prevailing winds and tropical storms guided by the 
strengthened blocking high pressure and nurtured by the un- 
usually warm late-Eemian tropical sea surface temperatures 
(Cortijo et al., 1999), which would favor more powerful trop- 
ical storms (Emanuel, 1987). This enhanced meridional tem- 
perature gradient - warmer tropics and cooler high latitudes 
- was enhanced by low obliquity of Earth’s spin axis in the 
late Eemian. 

3.9.2 21st century storms 

If GHGs continue to increase rapidly and ice melt grows, 
our simulations yield shutdown or major slowdown of the 
AMOC in the 21st century, implying an increase in severe 
weather. This is shown by zonal-mean temperature and eddy 
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kinetic energy changes in simulations of Sects. 3. 3-3. 6 with 
and without ice melt (Fig. 21). Without ice melt, surface 
warming is largest in the Arctic (Fig. 21, left), resulting in 
a decrease in lower tropospheric eddy energy. However, the 
surface cooling from ice melt increases surface and lower tro- 
pospheric temperature gradients, and in stark contrast to the 
case without ice melt, there is a large increase in midlatitude 
eddy energy throughout the midlatitude troposphere. The in- 
crease in zonal-mean midlatitude baroclinicity (Fig. 21) is 
in agreement with the localized. North Atlantic-centered in- 
creases in baroclinicity found in the higher-resolution simu- 
lations of Jackson et al. (2015) and Brayshaw et al. (2009). 

Increased baroclinicity produced by a stronger tempera- 
ture gradient provides energy for more severe weather events. 
Many of the most significant and devastating storms in east- 
ern North America and western Europe, popularly known 
as superstorms, have been winter cyclonic storms, though 
sometimes occurring in late fall or early spring, that gener- 
ate near-hurricane-force winds and often large amounts of 
snowfall (Chapter 11, Hansen, 2009). Continued warming of 
low-latitude oceans in coming decades will provide a larger 
water vapor repository that can strengthen such storms. If 
this tropical warming is combined with a cooler North At- 
lantic Ocean from AMOC slowdown and an increase in mid- 
latitude eddy energy (Fig. 21), we can anticipate more se- 
vere baroclinic storms. Increased high pressure due to cooler 
high-latitude ocean (Fig. 20) can make blocking situations 
more extreme, with a steeper pressure gradient between the 
storm’s low-pressure center and the blocking high, thus driv- 
ing stronger North Atlantic storms. 

Freshwater injection into the North Atlantic and Southern 
oceans increases sea level pressure at middle latitudes and 
decreases it at polar latitudes (Figs. 20, S22), but the impact 
is different in the North Atlantic than in the Southern Ocean. 
In the Southern Ocean the increased meridional temperature 
gradient increases the strength of westerlies in all seasons at 
all longitudes. In the North Atlantic Ocean the increase in 
sea level pressure in winter slows the westerlies (Fig. 20). 
Thus instead of a strong zonal wind that keeps cold polar air 
locked in the Arctic, there is a tendency for a less zonal flow 
and thus more cold air outbreaks to middle latitudes. 

4 Earth’s climate history 

Earth’s climate history is our richest source of information 
about climate processes. We first examine the Eemian or 
MIS 5e period, the last time Earth was as warm as today, 
because it is especially relevant to the issue of rapid sea level 
rise and storms when ice sheets existed only on Greenland 
and Antarctica. A fuller interpretation of late-Eemian climate 
events, as well as projection of climate change in the An- 
thropocene, requires understanding mechanisms involved in 
Earth’s millennial climate oscillations, which we discuss in 
the following subsection. 
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4.1 Eemian interglacial period (marine isotope 
substage MIS 5e) 

We first discuss Eemian sea level (Sect. 4.1.1), especially ev- 
idence for rapid sea level rise late in the Eemian to +6-9 m 
relative to today’s sea level, and then evidence for strong late- 
Eemian storms (Sect. 4.1.2). We provide in the Supplement 
more detailed geologic analysis of data on Eemian sea level, 
because the rapid late-Eemian sea level rise relates to our 
expectation of likely near-future events if rapid global warm- 
ing continues. In Sect. 4.1.3 we present evidence from ocean 
sediment cores for strong late-Eemian cooling in the North 
Atlantic associated with shutdown of the Atlantic meridional 
overturning circulation (AMOC), and in Sect. 4.1.4 we show 
that Earth orbital parameters in the late Eemian were consis- 
tent with cooling in the North Atlantic and global sea level 
rise from Antarctic ice sheet collapse. 

4.1.1 Eemian sea level 

Eemian sea level is of special interest because Eemian cli- 
mate was little warmer than today. Masson-Delmotte et 
al. (2013) conclude, based on multiple data and model 
sources, that peak Eemian temperature probably was only 
a few tenths of a degree warmer than today. Yet Eemian 
sea level reached heights several meters above today’s level 
(Land et al., 1967; Chen et al., 1991; Neumann and Hearty, 
1996; Hearty et al., 2007; Kopp et al., 2009; Dutton and Lam- 
beck, 2012; O’Leary et al., 2013; Dutton et al., 2015). 

Change of sea level within the Eemian period is partic- 
ularly relevant to concerns about ice sheet stability and the 
potential for rapid sea level rise. Hearty et al. (2007) used 
data from 15 sites around the world to construct an Eemian 
sea level curve that had sea level rising in the early Eemian 
to +2-3 m (“+” indicates above today’s sea level), possibly 
falling in the mid-Eemian to near today’s sea level, rapidly 
rising in the late Eemian to +6-9 m, and then plummeting as 
Earth moved from the Eemian into the 100000-year glacial 
period preceding the Holocene. Evidence from a variety of 
sources supports this interpretation, as discussed in the Sup- 
plement. 

The most comprehensive analyses of sea level and paleo- 
climate storms are obtained by combining information from 
different geologic sources, each with strengths and weak- 
nesses. Coral reefs, for example, allow absolute U/Th dating 
with age uncertainty as small as 1-2 ky, but inferred sea lev- 
els are highly uncertain because coral grows below sea level 
at variable depths as great as several meters. Carbonate plat- 
forms such as Bermuda and the Bahamas, in contrast, have 
few coral reefs for absolute dating, but the ability of carbon- 
ate sediments to cement rapidly preserves rock evidence of 
short-lived events such as rapid sea level rise and storms. 

The important conclusion, that sea level rose rapidly in the 
late Eemian by several meters, to +6-9 m, is supported by 
records preserved in both the limestone platforms and coral 
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Figure 21. Simulated zonal-mean atmospheric temperature (°C) and eddy kinetic energy (m 2 s 2 ) in 2078-2082 relative to 1880-1920 for 
A1B scenario and A1B plus 2.5 m ice melt in each hemisphere. 


reefs. Figure 6 of Hearty and Kindler (1995), for example, 
based on Bermuda and Bahamas geological data from marine 
and eolian limestone, reveals the rapid late-Eemian sea level 
rise and fall. Based on the limited size of the notches cut in 
Bahamian shore during the rapid late-Eemian level rise and 
crest, Neumann and Hearty (1996) inferred that this period 
was at most a few hundred years. Independently, Blanchon 
et al. (2009) used coral reef “back-stepping” on the Yucatan 
Peninsula, i.e., movement of coral reef building shoreward as 
sea level rises, to conclude that sea level in the late Eemian 
jumped 2-3 m within an “ecological” period, i.e., within sev- 
eral decades. 

Despite general consistency among these studies, consid- 
erable uncertainty remains about absolute Eemian sea level 
elevation and exact timing of end-Eemian events. Uncertain- 
ties include effects of local tectonics and glacio-isostatic ad- 
justment (GIA) of Earth’s crust. Models of GIA of Earth’s 
crust to ice sheet loading and unloading are increasingly used 
to improve assessments. O’Leary et al. (2013) use over 100 
corals from reefs at 28 sites along the 1400 km west coast of 
Australia, incorporating minor GIA corrections, to conclude 
that sea level in most of the Eemian was relatively stable at 
+3-4 m, followed by a rapid late-Eemian sea level rise to 
about +9 m. U-series dating of the corals has peak sea level 
at 118.1 ±1.4 kyb2k. 

A more complete discussion of data on Eemian sea level 
is provided in the Supplement. 

Late-Eemian sea level rise may seem a paradox, because 
orbital forcing then favored growth of Northern Hemisphere 
ice sheets. We will find evidence, however, that the sea level 
rise and increased storminess are consistent, and likely re- 
lated to events in the Southern Ocean. 

4.1.2 Evidence of end-Eemian storms in Bahamas and 
Bermuda 

Geologic data indicate that the rapid end-Eemian sea level 
oscillation was accompanied by increased temperature gradi- 
ents and storminess in the North Atlantic region. We summa- 


rize several interconnected lines of evidence for end-Eemian 
storminess, based on geological studies in Bermuda and the 
Bahamas referenced below. It is important to consider all the 
physical evidence of storminess rather than exclusively the 
transport mechanism of the boulders; indeed, it is essential 
to integrate data from obviously wave-produced runup and 
chevron deposits that exist within a few kilometers on North 
Eleuthera, Bahamas, as well as across the Bahamas Platform. 

The Bahama Banks are flat, low-lying carbonate plat- 
forms that are exposed as massive islands during glacials and 
largely inundated during interglacial high stands. From a tec- 
tonic perspective, the platforms are relatively stable, as indi- 
cated by near-horizontal ±2-3 m elevation of Eemian reef 
crests across the archipelago (Hearty and Neumann, 2001). 
During MIS 5e sea level high stands, an enormous volume of 
aragonitic oolitic grains blanketed the shallow, high-energy 
banks. Sea level shifts and storms formed shoals, ridges, 
and dunes. Oolitic sediments indurated rapidly (~ 10 1 to 
~ 10 2 years) once stabilized, preserving detailed and delicate 
lithic evidence of these brief, high-energy events. This shift- 
ing sedimentary substrate across the banks was inimical to 
coral growth, which partially explains the rarity of reefs dur- 
ing late MIS 5e. 

The preserved regional stratigraphic, sedimentary and ge- 
omorphic features attest to a turbulent end-Eemian transi- 
tion in the North Atlantic. As outlined below, a coastal gra- 
dient of sedimentological features corresponds with coastal 
morphology, distance from the coast, and increasing eleva- 
tion, reflecting the attenuating force and inland “reach” of 
large waves, riding on high late-Eemian sea levels. On rocky, 
steep coasts, giant limestone boulders were detached and cat- 
apulted onto and over the coastal ridge by ocean waves. On 
higher, Atlantic-facing built-up dune ridges, waves ran up to 
over 40 m elevation, leaving meter-thick sequences of fen- 
estral beds, pebble lenses, and scour structures. Across kilo- 
meters of low-lying tidal inlets and flats, “nested” chevron 
clusters were formed as stacked, multi-meter thick, tabular 
fenestral beds. 


www.atmos-chem-phys.net/16/3761/2016/ 


Atmos. Chem. Phys., 16, 3761-3812, 2016 


3782 





Figure 22. Megaboulders #1 (left) and #2 resting on MIS 5e eo- 
lianite at the crest of a 20 m high ridge with person (1.7 m) showing 
scale and orientation of bedding planes in the middle Pleistocene 
limestone. The greater age compared to underlying strata and dis- 
orientation of the primary bedding beyond natural in situ angles in- 
dicates that the boulders were wave-transported. 

The complexity of geomorphology and stratigraphy of 
these features are temporal measures of sustained sea level 
and storm events, encompassing perhaps hundreds of years. 
These features exclude a single wave cluster from a lo- 
cal point-source tsunami. Here we present data showing the 
connections among the megaboulders, runup deposits, and 
chevron ridges. 

Megaboulders 

In North Eleuthera enormous boulders were plucked from 
seaward middle Pleistocene outcrops and washed onto a 
younger Pleistocene landscape (Hearty and Neumann, 2001). 
The average lOOOt megaclasts provide a metric of powerful 
waves at the end of MIS 5e. Evidence of transport by waves 
includes that (1) they are composed of recrystallized oolitic- 
peloidal limestone of MIS 9 or 11 age (300-400 ky; Kindler 
and Hearty, 1996) and hammer-ringing hardness; (2) they 
rest on oolitic sediments typical of early to mid-MIS 5e that 
are soft and punky under hammer blows; (3) Cerion land 
snail fossils beneath boulder #4 (Hearty, 1997) correlate with 
the last interglacial period (Garrett and Gould, 1984; Hearty 
and Kaufman, 2009); (4) calibrated amino acid racemization 
(AAR) ratios (Hearty, 1997; Hearty et ah, 1998; Hearty and 
Kaufman, 2000, 2009) confirm the last interglacial age of the 
deposits as well as the stratigraphic reversal; (5) dips of bed- 
ding planes in boulders between 50 and 75° (Fig. 22) far ex- 
ceed natural angles; and (6) some of the largest boulders are 
located on MIS 5e deposits at the crest of the island’s ridge, 
proving that they are not karstic relicts of an ancient land- 
scape (Mylroie, 2008). 


J. Hansen et al.: Ice melt, sea level rise and superstorms 

The ability of storm waves to transport large boulders 
is demonstrated. Storms in the North Atlantic tossed boul- 
ders as large as 80 1 to a height +11 m on the shore on Ire- 
land’s Aran Islands (Cox et ah, 2012), this specific storm 
on 5 January 1991 being driven by a low-pressure system 
that recorded a minimum 946 mb, producing wind gusts to 
80 kn and sustained winds of 40 kn for 5 h (Cox et ah, 2012). 
Typhoon Haiyan (8 November 2013) in the Philippines pro- 
duced longshore transport of a 1 80 1 block and lifted boul- 
ders of up to ~ 24 1 to elevations as high as 10 m (May et ah, 
2015). May et ah (2015) conclude that these observed facts 
“demand a careful re-evaluation of storm-related transport 
where it, based on the boulder’s sheer size, has previously 
been ascribed to tsunamis”. 

The situation of the North Eleuthera megaboulders is spe- 
cial in two ways. First, all the large boulders are located at 
the apex of a horseshoe-shaped bay that would funnel en- 
ergy of storm waves coming from the northeast, the direction 
of prevailing winds. Second, the boulders are above a ver- 
tical cliff at right angles to the incoming waves, a situation 
that allows constructive interference of reflected and incom- 
ing waves (Cox et ah, 2012). The ability of waves hitting that 
cliff to produce large near-vertical splash is shown by a pho- 
tograph in the Supplement taken on 31 October 1991, when a 
storm in the North Atlantic produced large waves impacting 
Eleuthera. 

It is generally accepted that the boulders were wave- 
transported in the late Eemian. The boulders were deposited 
near complex chevron ridges and widespread runup deposits, 
which must be considered in analyzing wave-generating 
mechanisms. Lower elevation areas such as tidal inlets would 
have been flooded and scoured by the same waves, forming 
chevron ridges, and such large waves would also wash up 
onto higher, older ridges. 

Runup deposits 

Across several hundred kilometers of the Bahama Islands, 
older built-up dune ridges are mantled with wave runup de- 
posits that reach heights over +40 m (Fig. 23). They are gen- 
erally 1-5 m thick, fenestrae-filled, and seaward-sloping tab- 
ular beds (Wanless and Dravis, 1989; Chen et al., 1991; Neu- 
mann and Hearty, 1996; Tormey and Donovan, 2015). These 
stratigraphically youngest Eemian deposits mantle older MIS 
5e dune deposits on the shore-parallel ridges, and are the up- 
land correlative to wave-generated boulders and chevron for- 
mations. 

If these are deposits of powerful storms driven by an un- 
usually warm tropical ocean and strong temperature gradi- 
ents in the North Atlantic, as opposed to a tsunami, should 
there not be evidence of comparable end-Eemian storms in 
Bermuda? Indeed, along several kilometers of the north coast 
of Bermuda (Land et al., 1967; Vacher and Rowe, 1997; 
Hearty et al., 1998) there are seaward sloping planar beds 
rising to about +20 m. Although interpretations of these beds 
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vary, they are filled with beach fenestrae and stratigraphically 
of latest MIS 5e carbonate sediments equivalent to runup in 
the Bahamas. These planar beds contrast with older MIS 5e 
sedimentary (dune) structures that underlie them (Hearty et 
al., 1998). Massive subtidal cross beds comprise the seaward 
facies of the elevated beach beds, pointing to an exceptional 
energy anomaly on the normally tranquil, shallow, broad and 
protected north shore platform of Bermuda. 

Chevrons 

In the Bahama Islands, extensive oolitic sand ridges with a 
distinctive landward-pointing V shape are common, standing 
~ 5-15 m high across several kilometers on broad, low-lying 
platforms or ramps throughout the Atlantic-facing, deep- 
water margins of the Bahamas (Hearty et al., 1998). Hearty 
et al. (1998) examined 35 areas with chevron ridges across 
the Bahamas, which all point generally in a southwest direc- 
tion (S65° W) with no apparent relation to the variable aspect 
of the coastline, nor to a point source event as would be ex- 
pected for a tsunami generated by flank margin collapse. 

These chevron formations are the lowland correlative 
to the wave-generated rocky coast boulder deposits. The 
chevron ridges often occur in nested groups of several ridges 
(e.g.. North Eleuthera and Great Exuma; Hearty et al., 1998) 
and show multiple complex sets and subsets of fenestrae- 
filled beds, indicating the passage of a sustained interval of 
time late in the interglacial. Their definitive and complex 
characteristics preclude formation during a single tsunami 
event. 

The character of fenestral beds in both the Eemian chevron 
ridges and runup deposits change with increasing elevation 
and distance from shore, as do the abundance and geometry 
of fenestral pores (Tormey and Donovan, 2015): (1) at low 
elevations and in proximal locations, the chevron ridges are 
dominated by multiple truncated, thick, tabular, fenestrae- 
rich beds (Fig. 24a, b); (2) at moderate elevations and further 
inland, fenestrae are concentrated in discrete packages within 
eolianites, often associated with scour (Fig. 23) and rip-up 
clasts; and (3) in the highest and most distal eolian ridges, 
only rare, thin, discontinuous fenestrae beds can be found 
(Fig. 24c of Tormey and Donovan, 2015). This spatial tran- 
sition is improbable if torrential rain was falling across the 
area during a storm as asserted by Bain and Kindler (1994); 
rather this is exactly the pattern expected as waves attenuate 
with greater distance and elevation inland. 

Presence of a few eolian structures (Engel et al., 2015) 
does not imply that the chevron ridges are parabolic dunes; 
it suggests the deposits were sub-aerially exposed and wind 
blew during periods of relative quiescence (as commonly ob- 
served on today’s beaches after a storm). Unlike parabolic 
dunes, the chevron ridges are dominated by thick, low- 
angle (< 10°) seaward-dipping, aggradational oolitic bed- 
ding (Hearty et al., 1998; Tormey, 1999; Fig. 24a-c). Fore- 
set beds, diagnostic of migrating parabolic dunes, are rare 
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Figure 23. Photograph of runup deposits in a road cutting above 
+23 m ( 1 m scale in photo) on Old Land Road, Great Exuma Island, 
situated deep in Exuma Sound ~ 200 km south of North Eleuthera. 
The older built-up eolianite forms the lower half of the image; the 
upper half has multiple “packages” of planar, fenestrae-filled beach 
sets. The upper progression of sedimentary packages (labeled a- 
e) clearly shows an onlapping, rising sequence of beds, indicating 
increasing wave energy and degree of runup. Further, the individ- 
ual laminae of scour structures (arrows and inset image) display 
the same onlapping, upward-climbing succession. It would be im- 
possible to achieve such bedding if rain-saturated sediments were 
sloughing downhill on low-angle slopes under the influence of grav- 
ity, especially near the crest of a ridge. 


or absent from many chevron ridges, supporting formation 
primarily by waves (Hearty et al., 2002). Furthering the dis- 
tinction, fenestral porosity in low-angle bedding is prevalent 
throughout chevron ridges, occurring in repeated cycles of 
centimeter-thick beds that onlap the underlying strata, and 
often comprise meter-thick fenestrae-rich packages that can 
be followed in outcrop for tens of meters (Fig. 24a, b). 

Summary of evidence 

Alternative interpretations of the geologic data have been 
made (Bain and Kindler, 1994; Kindler and Strasser, 2000, 
2002; Engel et al., 2015); specifically, (1) the boulders were 
thrown by a tsunami caused by flank margin collapse in 
North Eleuthera, (2) beach fenestrae in runup and chevron 
beds were caused by heavy rainfall, and (3) the chevron 
beach ridges are parabolic dunes. These views are challenged 
by Hearty et al. (2002) and again here for the following rea- 
sons. (1) extensive research in the Bahamas has revealed no 
geologic evidence of a point-source tsunami radiating from 
North Eleuthera. A slow speed margin failure is possible, 
without a tsunami, and indeed such a flank margin collapse 
could have been initiated by massive storm waves impact- 
ing an over-steepened margin. (2) If heavy rainfall was a sig- 
nificant process in the formation of fenestrae in dunes, they 
should commonly occur in all dunes of all ages, which is not 
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Figure 24. Cross-section diagrams (Tormey, 1999) of Eemian chevron and dune deposits in North Eleuthera (a, b ~ 10 km west of megaboul- 
ders) showing geometry of bedding, fenestral porosity (lines of blue dots), and fossil roots (vertical wavy lines), (a) Chevron ridge exposure 
at Licrish Hill characterized by rising sequences of thick, tabular fenestral beds, (b) Chevron ridge exposure at Airport Junction characterized 
by rising sequences of thick, tabular fenestral beds, (c) Eolian ridge exposure at a higher elevation road cutting at Annie Bight (6 km south 
of megaboulders) characterized by dominantly backset and topset bedding with scattered, thin, wispy fenestrae beds. 


the case. (3) Carbonate dunes, particularly oolitic ones, gen- 
erally do not migrate unless exposed to extremely arid cli- 
mates, which contradicts point 2, and chevrons lack the most 
diagnostic feature of migration - foreset bedding. 

It is too random and chronologically coincidental to argue 
that the trilogy of evidence - boulders, runup deposits, and 
chevron ridges - was caused by unconnected processes. If 
large, long-period waves lifted 1000 1 boulders onto and over 
the coastal ridge, as is generally agreed, the same waves must 
have also impacted large areas of the eastern Bahamas, for 
which there is abundant documentation. A radiating pattern 
of landforms outward from a North Eleuthera point source, 
as from a tsunami generated from a local bank margin col- 
lapse, is not observed in the area or broader region. Absence 
of evidence for tsunamis on the United States East Coast re- 
futes the possibility of a large remote tsunami source. 

Our interpretation of these features is the most parsimo- 
nious, and we have argued that it is most consistent with 
the data. A common, synchronous, and non-random set of 
super-storm-related processes best explains boulder transport 
by waves, emplacement of runup deposits on older built-up 
ridges, and the formation of complex chevron deposits over 
time across lower areas of the Bahamas. Indeed, given the ge- 
ologic evidence of high seas and storminess from Bermuda 
and the Bahamas, Hearty and Neumann (2001) suggested 
“steeper pressure, temperature, and moisture gradients adja- 
cent to warm tropical waters could presumably spawn larger 
and more frequent cyclonic storms in the North Atlantic than 
those seen today”. 


We now seek evidence about end-Eemian climate change 
to help clarify how North Atlantic storms could have dis- 
persed such strong long-period, well-organized waves to the 
southwest. 

4.1.3 End-Eemian cold event: evidence from North 
Atlantic sediment cores 

We present here evidence for rapid end-Eemian cooling in 
the North Atlantic at a time with the tropics warmer than 
today. The cooling marked initial descent from interglacial 
conditions toward global ice age conditions, occurring at 
~ 118ky b2k in ocean cores with uncertainty ~2ky. It is 
identified by Chapman and Shackleton (1999) as cold event 
C26 in Greenland ice cores. 

This section discusses ocean core data, but we first note 
the relation with ice core data and dating uncertainties. Ice 
cores have become of great value for climate studies, partly 
because the relative timing of events in ice cores at differ- 
ent locations can be determined very accurately via marker 
events such as volcanic eruptions and CHq fluctuations, even 
though the absolute dating error in ice cores is comparable to 
the dating uncertainty in ocean cores. 

C26 is the cold phase of Dansgaard-Oeschger (D-O) cli- 
mate oscillation D-O 26 in the NGRIP (North Greenland 
Ice Core Project) ice core (NGRIP, 2004). C26 begins with 
a sharp cooling at 119.14ky b2k on the GICC05modelext 
timescale (Rasmussen et al., 2014). The GICC05 timescale 
is based on annual layer counting in Greenland ice cores 
for the last 60 ky and an ice-flow-model extension for ear- 
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Ocean and ice core sites; Surface and deep currents 



Figure 25. Ocean and ice core sites and simplified sketch of up- 
per ocean North Atlantic Current and North Atlantic Deep Water 
return flow. In interglacial periods the North Atlantic Current ex- 
tends further north, allowing the Greenland-Iceland-Norwegian Sea 
to become an important source of deepwater formation. 


Her times. An alternative timescale is provided by Antarc- 
tic ice core chronology AICC2012 (Bazin et al., 2013; 
Veres et al., 2013), on which Greenland ice core records 
are synchronized via global markers, mainly oscillations of 
atmospheric CH4 amount, which is globally well mixed. 
C26 on Greenland is at 1 16.72 ky b2k on the AICC2012 
timescale. Figure S19 shows the difference between GICC05 
and AICC2012 timescales versus time. 

This age uncertainty for C26 is consistent with the ice 
core 2a error estimate of 3.2 ky at Eemian time (Bazin et 
al., 2013). Despite this absolute age uncertainty, we can use 
Greenland data synchronized to the AICC2012 timescale to 
determine the relative timing of Greenland and Antarctic cli- 
mate changes (Sect. 4.2.1) to an accuracy of a few decades 
(Bazin et al., 2013). 

Sediment cores from multiple locations provide informa- 
tion not only on ocean temperature and circulation (Fig. 25) 
but also on ice sheet via information inferred from ice-rafted 
debris. Comparison of data from different sites is affected by 
inaccuracy in absolute dating and use of different age mod- 
els. Dating of sediments is usually based on tuning to the 
timescale of Earth orbital variations (Martinson et al., 1987) 
or “wiggle matching” to another record (Sirocko et al., 2005), 
which limits accuracy to several thousand years. Temporal 
resolution is limited by bioturbation of sediments; thus res- 
olution varies with core location and climate (Keigwin and 
Jones, 1994). For example, high deposition rates during ice 
ages at the Bermuda Rise yield a resolution of a few decades, 
but low sedimentation rates during the Eemian yield a reso- 
lution of a few centuries (Lehman et al., 2002). Lateral trans- 
port of sedimentary material prior to deposition complicates 
data interpretation and can introduce uncertainty, as argued 
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specifically regarding data from the Bermuda Rise (Ohk- 
ouchi et al., 2002; Engelbrecht and Sachs, 2005). 

Adkins et al. (1997) analyzed a sediment core (MD95- 
2036; 34° N, 58° W) from the Bermuda Rise using an age 
model based on Martinson et al. (1987) orbital tuning with 
the MIS stage 5-6 transition set at 131 ky b2k and the stage 
5d-5e transition at 1 14 ky b2k. They found that oxygen iso- 
tope <5 18 0 of planktonic (near-surface dwelling) foraminifera 
and benthic (deep-ocean) foraminifera both attain full inter- 
glacial values at ~ 128 ky b2k and remain nearly constant 
for ~ 10 ky (their Fig. 2). Adkins et al. (1997) infer that “late 
within isotope stage 5e (~ 118 ky b2k), there is a rapid shift 
in oceanic conditions in the western North Atlantic. . . ”, They 
And in the sediments at that point an abrupt increase in clays 
indicative of enhanced land-based glacier melt and an in- 
crease in high nutrient “southern source waters”. The latter 
change implies a shutdown or diminution of NADW forma- 
tion that allows Antarctic Bottom Water (AABW) to push 
into the deep North Atlantic Ocean (Duplessy et al., 1988; 
Govin et al., 2009). Adkins et al. (1997) continue: “The rapid 
deep and surface hydrographic changes found in this core 
mark the end of the peak interglacial and the beginning of 
climate deterioration towards the semi-glacial stage 5d. Be- 
fore and immediately after this event, signaling the impeding 
end of stage 5e, deep-water chemistry is similar to modern 
NADW.” This last sentence refers to a temporary rebound to 
near-interglacial conditions. In Sect. 4.2.4 we use accurately 
synchronized Greenland and Antarctic ice cores, which also 
reveal this temporary end-Eemian climate rebound, to inter- 
pret the glacial inception and its relation to ice melt and late- 
Eemian sea level rise. 

Ice-rafted debris (IRD) found in ocean cores provides a 
useful climate diagnostic tool (Heinrich, 1988; Hemming, 
2004). Massive ice rafting (“Heinrich”) events are often as- 
sociated with decreased NADW production and shutdown or 
slowdown of the Atlantic meridional overturning circulation 
(AMOC) (Broecker, 2002; Barreiro et al., 2008; Srokosz et 
al., 2012). However, ice rafting occurs on a continuum of 
scales, and significant IRD is found in the cold phase of all 
the 24 D-O climate oscillations first identified in Greenland 
ice cores (Dansgaard et al., 1993). D-0 events exhibit rapid 
warming on Greenland of at least several degrees within a 
few decades or less, followed by cooling over a longer pe- 
riod. Chapman and Shackle ton (1999) found IRD events in 
the NEAP18K core for all D-O events (C19-C24) within the 
core interval that they studied, and they also labeled two ad- 
ditional events (C25 and C26). C26 did not produce identifi- 
able IRD at the NEAP18K site, but it was added to the series 
because of its strong surface cooling. 

Lehman et al. (2002) quantify the C26 cooling event us- 
ing the same Bermuda Rise core (MD95-2036) and age 
model as Adkins et al. (1997). Based on the alkenone paleo- 
temperature technique (Sachs and Lehman, 1999), Lehman 
et al. (2002) find a sharp SST decrease of ~ 3 °C (their Fig. 1) 
at ~ 1 18 ky BP, coinciding with the end-Eemian shoulder of 
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the benthic <$ 18 0 plateau that defines stage 5e in the deep 
ocean. The SST partially recovered after several centuries, 
but C26 marked the start of a long slide into the depths of 
stage 5d cold, as ice sheets grew and sea level fell ~ 50 m in 
10 ky (Lambeck and Chappell, 2001; Rohling et ah, 2009). 
Lehman et al. (2002) wiggle-match the MD95-2036 and 
NEAP18K cores, finding a simple adjustment to the age 
model of Chapman and Shackleton (1999) that maximizes 
correlation of the benthic <5 ls O records with the Adkins et 
al. (1997) (5 ls O record. Specifically, they adjust the NEAP 
timescale by +4ky before the MIS 5b /> 18 0 minimum and 
by +2 ky after it, which places C26 cooling at 1 18 ky b2k in 
both records. They give preference to the Adkins et al. (1997) 
age scale because it employs a 230 Th-based timescale be- 
tween 100 and 130 ky b2k. 

We do not assert that the end-Eemian C-26 cooling was 
necessarily at 1 18 ky b2k, but we suggest that the strong rapid 
cooling observed in several sediment cores in this region of 
the subtropical and midlatitude North Atlantic Drift at about 
this time were all probably the same event. Such a large cool- 
ing lasting for centuries would not likely be confined to a 
small region. The dating models in several other studies place 
the date of the end-Eemian shoulder of the deep-ocean 5 ls O 
and an accompanying surface cooling event in the range 1 16— 
118kyb2k. 

Kandiano et al. (2004) and Bauch and Kandiano (2007) 
analyze core M23414 (53° N, 17° W), west of Ireland, find- 
ing a major SST end-Eemian cooling that they identify as 
C26 and place at 117ky b2k. The Iky change in the tim- 
ing of this event compared with Lehman et al. (2002) is 
due to a minor change in the age model; specifically, Bauch 
and Kandiano say that “the original age model of MD95- 
2036 (Lehman et al., 2002) has been adjusted to our core 
M23414 by alignment of the 4 per mil level in the benthic 
(5 ls O records (at 130 ka in M23414) and the prominent C24 
event in both cores”. Bauch and Erlenkeuser (2008) and H. 
Bauch et al. (2012) examine ocean cores along the North At- 
lantic Current including its continuation into the Nordic Seas. 
They find that, in the Greenland-Iceland-Norwegian (GIN) 
seas, unlike middle latitudes, the Eemian was warmest near 
the end of the interglacial period. The age model employed 
by Bauch and Erlenkeuser (2008) has the Eemian about 2ky 
younger than the Adkins et al. (1997) age model, Bauch 
and Erlenkeuser (2008) having the benthic <5 ls O plateau at 
~1 16-124 ky BP (their Fig. 6). Rapid cooling they illus- 
trate there at ~ 116.6kyBP for core M23071 on the Vpring 
Plateau (67° N, 3° E) likely corresponds to the C26 end- 
Eemian cooling event. 

Identification of end-Eemian cooling in ocean cores is 
hampered by the fact that Eemian North Atlantic climate 
was more variable than in the Holocene (Fronval and Jansen, 
1996). There were at least three cooling events within the 
Eemian, each with minor increases in IRD, which are labeled 
C27, C27a and C27b by Oppo et al. (2006); see their Fig. 2 
for core site ODP-980 in the eastern North Atlantic (55° N, 


J. Hansen et al.: Ice melt, sea level rise and superstorms 

15° W) near Ireland. High- (sub-centennial) resolution cores 
in the Eirik drift region (MD03-2664, 57° N, 49° W) near the 
southern tip of Greenland reveal an event with rapid cooling 
accompanied by reduction in NADW production (Irvali et al., 
2012; Galaasen et al., 2014), which they place at ~ 117ky 
b2k. However, their age scale has the benthic 5 18 0 shoulder 
at ~ 115ky b2k (Fig. SI of Galaasen et al., 2014), so that 
event may have been C27b, with C26 being stronger cooling 
that occurred thereafter. 

4.1.4 Eemian timing consistency with insolation 
anomalies 

Glacial-interglacial climate cycles are related to insolation 
change, as shown persuasively by Hays et al. (1976). Each 
“termination” (Broecker, 1984) of glacial conditions in the 
past several hundred thousand years coincided with a large 
positive warm-season insolation anomaly at the latitude of 
North American and Eurasian ice sheets (Raymo, 1997; Pail- 
lard, 2001). The explanation is that positive summer inso- 
lation anomalies (negative in winter) favor increased sum- 
mer melting and reduced winter snowfall, thus shrinking ice 
sheets. 

Termination timing is predicted better by high Northern 
Hemisphere late spring (April-May-June) insolation than by 
summer anomalies. For example, Raymo (1997) places mid- 
points of Termination I and II (preceding the Holocene and 
Eemian) at 13.5 and 128-131 ky b2k. Late spring insola- 
tion maxima are at 13.2 and 129.5 ky b2k (Fig. 26a). The 
AICC2012 ice core chronology (Bazin et al., 2013) places 
Termination II at 128.5 ky b2k, with 2er uncertainty of 3.2 ky. 
Late spring irradiance maximizes warm-season ice melt by 
producing the earliest feasible warm-season ice sheet dark- 
ening via snow melt and snow recrystallization (Hansen et 
al., 2007b). 

Summer insolation anomalies are also shown in Fig. 26, 
because interglacial periods can be expected to continue as 
long as summer insolation is large enough to prevent ice 
sheet genesis. Summer insolation anomalies at 60° N became 
negative at ~118ky b2k (Fig. 26a). Dating of insolation 
anomalies has high absolute accuracy, unlike ocean and ice 
core dating, as orbital anomalies are based on well-known 
planetary orbital mechanics (Berger, 1978). 

Late Eemian sea level rise might appear to be a paradox, 
because glacial-interglacial sea level change is mainly a re- 
sult of the growth and decay of Northern Hemisphere ice 
sheets. Northern warm-season insolation anomalies were de- 
clining rapidly in the latter part of the Eemian (Fig. 26a), so 
Northern Hemisphere ice should have been just beginning to 
grow. We suggest that the explanation for a late-Eemian sea 
level maximum is a late-Eemian collapse of Antarctic ice fa- 
cilitated by the positive warm-season insolation anomaly on 
Antarctica and the Southern Ocean during the late Eemian 
(Fig. 26b) and possibly aided by an AMOC shutdown, which 
would increase warming of the Southern Ocean. 


Atmos. Chem. Phys., 16, 3761-3812, 2016 


www.atmos-chem-phys.net/16/3761/2016/ 



J. Hansen et al.: Ice melt, sea level rise and superstorms 


3787 


(a) Insolation anomalies at 60° N (b) Insolation anomalies at 60° S 



Figure 26. Summer (June-July-August) and late spring (April-May-June) insolation anomalies (W m -) at 60° N and summer (December- 
January-February) and late spring (October-November-December) anomalies at 60° S. 


Persuasive evidence for this interpretation is provided by 
detailed paleoclimate data discussed in the next section, and 
is supported by modeling of relevant climate mechanisms. 
We will show that these mechanisms in turn help to explain 
ongoing climate change today, with implications for contin- 
uing climate change this century. 

4.2 Millennial climate oscillations 

Paleoclimate data are essential for understanding the ma- 
jor climate feedbacks. Processes of special importance are 

(1) the role of the Southern Ocean in ventilating the deep 
ocean, affecting CO 2 control of global temperature, and 

(2) the role of subsurface ocean warming in ice shelf melt, 
affecting ice sheet disintegration and sea level rise. An under- 
standing of timescales imparted by the ocean and the carbon 
cycle onto climate change is important, so that slow paleo- 
ice-sheet changes are not ascribed to ice physics, when the 
timescale is actually set elsewhere. 

Major glacial-interglacial climate oscillations are spurred 
by periodic variation of seasonal and geographical insola- 
tion (Hays et al., 1976). Insolation anomalies are due to 
slow changes of the eccentricity of Earth’s orbit, the tilt 
of Earth’s spin axis, and the precession of the equinoxes, 
and thus the day of year at which Earth is closest to the 
Sun, with dominant periodicities near 100000, 40000 and 
20000 years (Berger, 1978). These periods emerge in long 
climate records, yet a large fraction of climate variability at 
any site is stochastic (Wunsch, 2004; Lisiecki and Raymo, 
2005). Such behavior is expected for a weakly-forced system 
characterized by amplifying feedbacks, complex dynamics, 
and multiple sources of inertia with a range of timescales. 

Large glacial-interglacial climate change and stochastic 
variability are a result of two strong amplifying feedbacks, 
surface albedo and atmospheric CO 2 . Orbit-induced insola- 
tion anomalies, per se, cause a direct climate forcing, i.e., 
an imposed Earth energy imbalance, only of the order of 
0.1 W m -2 , but the persistent regional insolation anomalies 
spur changes of ice sheet size and GHGs. The albedo and 


GHG changes arise as slow climate feedbacks, but they are 
the forcings that maintain a quasi-equilibrium climate state 
nearly in global radiative balance. 

Glacial-interglacial albedo and greenhouse forcings are 
each ~3W m -2 (Fig. 27e, f) 2 . These forcings fully account 
for glacial-interglacial global temperature change with a cli- 
mate sensitivity of 0.5-1 °C per W m -2 (Hansen et al., 2008; 
Masson-Delmotte et al., 2010; Palaeosens, 2012). 

The insolation anomaly peaking at 129.5 ky b2k (Fig. 27a) 
succeeded in removing ice sheets from North America and 
Eurasia and in driving atmospheric CO 2 up to ~ 285 ppm, 
as discussed below. However, smaller climate oscillations 
within the last glacial cycle are also instructive about ice 
feedbacks. Some of these oscillations are related to weak 
insolation anomalies, and all are affected by predominately 
amplifying climate feedbacks. 

Insolation anomalies peaking at 107 and 86 ky b2k 
(Fig. 27a) led to ~ 40 m sea level rises at rates exceeding 
lmcentury -1 (Stirling et al., 1998; Cutler et al., 2003) in 
early MIS 5c and 5a (Fig. 27f), but CO 2 did not rise above 
250 ppm and interglacial status (with large ice sheets only 
on Greenland and Antarctica) was not achieved. CO 2 then 
continued on a 100 ky decline until ~ 18 ky b2k. Sea level 
continued its long decline, in concert with CO 2 , reaching a 
minimum at least 120 m below today’s sea level (Peltier and 
Fairbanks, 2006; Lambeck et al., 2014). 

Progress achieved by the paleoclimate and oceanographic 
research communities allows interpretation of the role of the 
Southern Ocean in the tight relationship between CO 2 and 


2 Other parts of Fig. 27 are discussed later, but they are most in- 
formative if aligned together. In interpreting Fig. 27, note that long- 
lived greenhouse gas amounts in ice cores have global relevance, but 
ice core temperatures are local to Greenland and Antarctica. Also, 
because our analysis does not depend on absolute temperature, we 
do not need to convert the temperature proxy, 5 18 0, into an esti- 
mated absolute temperature. We include CH 4 and N 2 O in the total 
GHG climate forcing, but we do not discuss the reasons for CH 4 
and N 2 O variability (see Schilt et al., 2010), because CO 2 provides 
~ 80 % of the GHG forcing. 
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(a) Late spring insolation anomaly at 60° N and 60° S 




(f) Sea level and albedo forcing 



Figure 27. (a) Late spring insolation anomalies relative to the mean for the past million years, (b) S 18 Oi ce of composite Greenland ice cores 
(Rasmussen et al., 2014) with Heinrich events of Guillevic et al. (2014). (c, d) <5 1 8 Oi ce of EDML Antarctic ice core (Ruth et al., 2007), 
multi-ice-core CO 2 , CH 4 , and N 2 O based on a spline fit with a 1000-year cut-off (Schilt et al., 2010); scales are such that CO 2 and S 18 0 
means coincide and standard deviations have the same magnitude, (e) GHG forcings from equations in Table 1 of Hansen et al. (2000), but 
with the CO 2 , CH 4 , and N 2 O forcings multiplied by factors 1.024, 1.60, and 1.074, respectively, to account for each forcing’s “efficacy” 
(Hansen et al., 2005a), with CH 4 including a factor of 1.4 to account for indirect effect on ozone and stratospheric water vapor, (f) Sea 
level data from Grant et al. (2012) and Lambeck et al. (2014) and ice sheet model results from de Boer et al. (2010). Marine isotope stage 
boundaries from Lisiecki and Raymo (2005). Panels (b-e) are on AICC2012 timescale (Bazin et al.. 2013), 


temperature, as well as discussion of the role of subsurface 
ocean warming in sea level rise. Both topics are needed to 
interpret end-Eemian climate change and ongoing climate 
change. 

4.2.1 Southern Ocean and atmospheric CO 2 

There is ample evidence that reduced atmospheric CO 2 in 
glacial times, at least in substantial part, results from in- 
creased stratification of the Southern Ocean that reduces ven- 
tilation of the deep ocean (Toggweiler, 1999; Anderson et 


al., 2009; Skinner et al., 2010; Tschumi et al., 2011; Burke 
and Robinson, 2012; Schmitt et al., 2012; Marcott et al., 
2014). Today the average “age” of deep water, i.e., the time 
since it left the ocean surface, is ~ 1000 years (DeVries and 
Primeau, 2011), but it was more than twice that old during 
the Last Glacial Maximum (Skinner et al., 2010). The South- 
ern Ocean dominates exchange between the deep ocean and 
atmosphere because ~ 80 % of deep water resurfaces in the 
Southern Ocean (Lumpkin and Speer, 2007), as westerly cir- 
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cumpolar winds and surface flow draw up deep water (Talley, 

2013) . 

Mechanisms causing more rapid deep-ocean ventilation 
during interglacials include warmer Antarctic climate that in- 
creases heat flux into the ocean and buoyancy mixing that 
supports upwelling (Watson and Garabato, 2006), poleward 
shift of the westerlies (Toggweiler et al., 2006), and reduced 
sea ice (Keeling and Stephens, 2001). Fischer et al. (2010) 
question whether the latitudinal shift of westerlies is an im- 
portant contributor; however, the basic point is the empiri- 
cal fact that a warmer interglacial Southern Ocean produces 
faster ventilation of the deep ocean via a combination of 
mechanisms. 

Although a complete quantitative understanding is lack- 
ing for mechanisms to produce the large glacial-interglacial 
swings of atmospheric CO 2 , we can safely assume that deep- 
ocean ventilation acts to oppose the sequestration of carbon 
in the ocean by the “pumps” that move carbon from the 
surface to ocean depths, and changes in the ventilation rate 
have a significant effect on atmospheric CCK Ridgwell and 
Arndt (2015) describe several conceptual pumps: (1) the or- 
ganic carbon pump, in which the sinking material also con- 
trols nutrient cycling by the ocean; (2) the carbonate pump, 
with biological precipitation of mainly calcium carbonate, a 
fraction of which escapes dissolution to form a new geologi- 
cal carbon reservoir; (3) the simple solubility pump, as CO 2 
is more soluble in the cold polar waters where deep water 
forms; and (4) a microbial carbon pump that seems capable 
of altering deep-ocean dissolved organic carbon. 

No doubt the terrestrial biosphere also contributes to 
glacial-interglacial atmospheric CO 2 change (Archer et al., 
2000; Sigman and Boyle, 2000; Kohler et al., 2005; Men- 
viel et al., 2012; Fischer et al., 2015). Also, the efficacy 
of the ocean pumps depends on terrestrial conditions, e.g., 
dust-borne iron fertilization of the biological pump (Mar- 
tin and Fitzwater, 1988) contributes to millennial and full 
glacial CO 2 drawdown (Martinez-Garcia et al., 2014). More- 
over, the Southern Ocean is not the only conduit to the deep 
ocean; for example, AMOC changes are associated with at 
least two rapid CO 2 increases of about 10 ppm, as revealed 
by a high-resolution West Antarctic ice core (Marcott et al., 

2014) . Nevertheless, it is reasonable to hypothesize that se- 
questration of CO 2 in the glacial ocean is the largest cause of 
glacial-interglacial CO 2 change, and it is known that ocean 
ventilation occurs mainly via the Southern Ocean. 

Southern Ocean ventilation, as the dominant cause of 
atmospheric CO 2 change, helps explain temperature-C02 
leads and lags. Temperature and CO 2 rises are almost con- 
gruent at ice age terminations (Masson-Delmotte et al., 2010; 
Pedro et al., 2012; Parrenin et al., 2013). Southern Ocean 
temperature is expected to lead, spurring deep-ocean ven- 
tilation and atmospheric CO 2 increase, with global temper- 
ature following. Termination I is dated best and Shakun et 
al. (2012) have reconstructed global temperature for that pe- 
riod, finding evidence for this expected order of events. 
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Correlation of S ls O and CO 2 over the past 140 ky 
(Fig. 27c) is 84.4 %, with CO 2 lagging by 760 years. For the 
period 1 00-20 ky b2k, which excludes the two terminations, 
the correlation is 77.5%, with CO 2 lagging by 1040 years. 
Briefer lag for the longer period and longer lag during glacial 
inception are consistent with the rapid deep-ocean ventilation 
that occurs at terminations. 

4.2.2 CO 2 as a climate control knob 

CO 2 is the principal determinant of Earth’s climate state, the 
radiative “control knob” that sets global mean temperature 
(Lacis et al., 2010, 2013). Degree of control is shown by 
comparison of CO 2 amount with Antarctic temperature for 
the past 800000 years (Fig. 28a). Control should be even 
tighter for global temperature than for Antarctic tempera- 
ture, because of regional anomalies such as Antarctic tem- 
perature overshoot at terminations (Masson-Delmotte et al., 
2006, 2010), but global data are not available. 3 

The CO 2 dial must be turned to ~ 260 ppm to achieve a 
Holocene-level interglacial. CO 2 ~250ppm was sufficient 
for quasi-interglacials in the period 800-450 ky b2k, with 
sea level 10-25 m lower than in the Holocene (Fig. S18 of 
Hansen et al., 2008). Interglacials with CO 2 ~ 280ppm, i.e., 
the Eemian and Holsteinian (~400ky b2k), were warmer 
than the Holocene and had sea level at least several meters 
higher than today. 

CO 2 and albedo change are closely congment over the last 
800000 years (Fig. S18 of Hansen et al., 2008). GHG and 
albedo forcings, which are both amplifying feedbacks that 
boost each other, are each of amplitude ~3W m -2 . So why 
do we say that CO 2 is the control knob? 

First, CO 2 , in addition to being a slow climate feedback, 
changes independently of climate. Natural CO 2 change in- 
cludes an increase to ~ 1000 ppm about 50 million years ago 
(Zachos et al., 2001) as a result of plate tectonics, specifically 
volcanic emissions associated with movement of the Indian 
plate across the Tethys Ocean and collision with Asia (Kent 
and Muttoni, 2008). Humankind, mainly by burning fossil 
fuels, also moves the CO 2 control knob. 


3 The tight fit of CO 2 and Antarctic temperature (Fig. 28a) im- 
plies an equilibrium Antarctic sensitivity of 20 °C for 2 x CO 2 
(4Wnf 2 ) forcing (200 -*300 ppm forcing is ~2.3Wm -2 ; Ta- 
ble 1 of Hansen et al., 2000) and thus 10 °C global climate sensi- 
tivity (Antarctic temperature change is around twice that of global 
change) with CO 2 taken as the ultimate control knob, i.e., if 
snow/ice area and other GHGs are taken to be slaves to CO 2 -driven 
climate change. This implies a conventional climate sensitivity of 
4°C for 2 x CO 2 , as GHG and albedo forcings are similar for 
glacial-to-interglacial climate change and non-CCB GHGs account 
for ~ 20 % of the GHG forcing. The inferred sensitivity is reduced 
to 2.5-3 °C for 2 x CO 2 if, as some studies suggest, global mean 
glacial-interglacial temperature change is only about one-third of 
the Antarctic temperature change (Palaeosens, 2012; Hansen et al., 
2013b). 
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(a) Antarctic temperature change and C0 2 



Time (EDC3, ky b2k) 

Figure 28. (a) Antarctic (Dome C) temperature relative to last lOky (Jouzel et al., 2007) on AICC2012 timescale and CCD amount (Luthi 
et al., 2008). Temperature scale is such that standard deviation of T and CO 2 are equal, yielding AT (°C) = 0.114 ACO 2 (ppm), (b) Late 
spring insolation anomalies at 60° N and 60° S. 


Second, CCD is more recalcitrant than snow and ice, i.e., 
its response time is longer. CO 2 inserted into the climate sys- 
tem, by humans or plate tectonics, remains in the climate 
system of the order of 100000 years before full removal by 
weathering (Archer, 2005). Even CO 2 exchange between the 
atmosphere (where it affects climate) and ocean has a lag of 
the order of a millennium (Fig. 27). In contrast, correlations 
of paleo-temperatures and sea level show that lag of sea level 
change behind temperature is of the order of a century, not a 
millennium (Grant et al., 2012). 

We suggest that limitations on the speed of ice volume 
(and thus sea level) changes in the paleo-record are more a 
consequence of the pace of orbital changes and CCD changes, 
as opposed to being a result of lethargic ice physics. “Fast” 
changes of CO 2 have been identified, e.g., an increase of 
~ 10 ppm in about a century at ~ 39.6 ky b2k (Ahn et al., 
2012) and three increases of 10-15 ppm each within 1-2 cen- 
turies during the deglaciation following the last ice age (Mar- 
cott et al., 2014), but the magnitude of these CO 2 increases 
is not sufficient to provide a good empirical test of ice sheet 
sensitivity to the CCD forcing. 

Dominance of SMOC, the Southern Ocean meridional 
overturning circulation, in affecting the CO 2 control knob 
and thus glacial-interglacial change is contrary to the idea 
that the AMOC is a prime driver that flips global climate be- 
tween quasi-stable glacial and interglacial states, yet AMOC 
retains a significant role. AMOC can affect CO 2 via the vol- 
ume and residence time of NADW, but its largest effect is 
probably via its impact on the Southern Ocean. When AMOC 
is not shut down it cools the Southern Hemisphere, trans- 
ferring heat from the Southern to the Northern Hemisphere 
at a rate of ~ 1 PW, which is ~ 4 W m -2 averaged over a 
hemisphere (Crowley, 1992). However, the Southern Ocean 
slowly warms when AMOC shuts (or slows) down; the re- 
sponse time is of the order of 1000 years because of the 


Southern Ocean’s large thermal inertia (Stocker and Johnson, 
2003). These mechanisms largely account for the nature of 
the “bipolar seesaw” (Broecker, 1998; Stocker, 1998; Stenni 
et al., 2011; Landais et al., 2015), including the lag between 
AMOC slowdown and Antarctic warming. 

4.2.3 Dansgaard-Oeschger (D-O) events and 
subsurface ocean warming 

The magnitude and rapidity of Greenland climate change 
during D-0 events would deter prediction of human-made 
climate effects if D-O events remained a mystery. Instead, 
however, enough is now understood about D-O events that 
they provide insight related to the vulnerability of ice shelves 
and ice sheets, including the role of subsurface ocean warm- 
ing. 

Broecker (2000) inferred from the rapidity of D-O warm- 
ings that a reduction of sea ice cover was probably involved. 
Li et al. (2005, 2010) modeling showed that removal of 
Nordic Seas ice cover is needed to yield the magnitude of 
observed Greenland warming. The spatial gradient of D- 
O warming, with smaller warming in northwest Greenland, 
agrees with that picture (Guillevic et al., 2013; Buizert et 
al., 2014). Such sea ice change is consistent with changes in 
deuterium excess in Greenland ice cores at D-O transitions, 
which indicate shifts of Greenland moisture source regions 
(Masson-Delmotte et al., 2005; Jouzel et al., 2007). 

Fluckiger et al. (2006), Alvarez-Solas et al. (2010, 2011, 
2013) and Marcott et al. (2011) noted modern and paleo-data 
that point to ocean-ice shelf interaction as key to the ice 
discharge of accompanying Heinrich events, and they used 
a range of models to support this interpretation and over- 
turn earlier suggestions of a central role for ice sheets via 
binge-purge oscillations (MacAyeal, 1993) or outburst flood- 
ing from subglacial reservoirs (Alley et al., 2006). Shaffer et 
al. (2004) and Petersen et al. (2013) conclude that subsur- 
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face ocean warming in the North Atlantic takes place dur- 
ing the stadial (cold) phase of all D-0 events, and eventually 
this subsurface warming leads to ice shelf collapse or retreat, 
ice rafting, sea level rise, and sea ice changes. Rasmussen et 
al. (2003) examined ocean cores from the southeast Labrador 
Sea and found that for all 1 1 Heinrich events “the icy surface 
water was overlying a relatively warm, poorly ventilated and 
nutrient rich intermediate water mass to a water depth of at 
least 1251 m”. Collapse of a Greenland ice shelf fronting the 
Jakoshavn ice stream during the Younger Dryas cold event 
has been documented (Rinterknecht et al., 2014), apparently 
due to subsurface warming beneath the ice shelf leading to 
rapid discharge of icebergs. 

Some D-0 details are uncertain, e.g., the relation between 
changing sea ice cover and changing location of deep wa- 
ter formation (Rahmstorf, 1994) and whether an ice shelf 
between Greenland and Iceland contributed to the sea ice 
variability (Petersen et al., 2013). However, ocean-ice inter- 
actions emerge as key mechanisms, spurred by subsurface 
ocean warming, as ocean stratification slows but does not 
stop northward heat transport by AMOC. 

We consider a specific D-O event for the sake of discussing 
mechanisms. D-0 22 cold phase, labeled C22 in ocean cores 
and coinciding with Heinrich H8 (Fig. 27), occurred as 
Northern Hemisphere insolation was rising (Fig. 27a). The 
North Atlantic surface was cooled by rapid ice discharge; 
sea level rose more than 40 m, a rate exceeding 1 .6 m per 
century (Cutler et al., 2003). Ice discharge kept the North 
Atlantic highly stratified, slowing AMOC. Antarctic warm- 
ing from a slowed AMOC increases almost linearly with the 
length of the D-0 cold phase (Fig. 3 of EPICA Community 
Members, 2006; Fig. 6 of Capron et al., 2010) because of 
the Southern Ocean’s large heat capacity (Stocker and John- 
son, 2003). Antarctic warming, aided by the 2500-year D-0 
22 event, spurred SMOC enough to raise atmospheric CO 2 
40 ppm (Fig. 27c). 

As the Antarctic warmed, ocean heat transport to the North 
Atlantic would have increased, with most heat carried at 
depths below the surface layer. When the North Atlantic be- 
came warm enough at depth, stratification of cold fresh sur- 
face water eventually could not be maintained. The warming 
breakthrough may have included change in NADW forma- 
tion location (Rahmstorf, 1994) or just large movement of 
the polar front. Surface warming east of Greenland removed 
most sea ice and Greenland warmed ~ 10 °C (Capron et al., 
2010). As the warm phase of D-O 21 began, AMOC was 
pumping heat from the Antarctic into the Nordic Seas and 
Earth must have been slightly out of energy balance, cooling 
to space, so both Antarctica and Greenland slowly cooled. 
Once the North Atlantic had cooled enough, sea ice formed 
east of Greenland again, ice sheets and ice shelves grew, sea 
level fell, and the polar front moved southward. 

Sea level rise associated with D-O events covers a wide 
range. Sea level increases as large as ~ 40 m were associated 
with large insolation forcings at 107 and 86 ky b2k (Fig. 27). 
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However, rapid sea level change occurred even when forc- 
ing was weak. Roche et al. (2004) conclude from analyses of 
<5 ls O that H4, at a time of little insolation forcing (~40ky 
b2k. Fig. 27), produced 1.9 ± 1.1 m sea level rise over 250 
± 150 years. Sea level rise as great as 10-15 m occurred in 
conjunction with some other D-0 events during 65-30 ky 
b2k (Lambeck and Chappell, 2001; Yokoyama et al., 2001; 
Chappell, 2002). 

Questions about possible D-O periodicity and external 
forcing were raised by a seeming 1470-year periodicity 
(Schulz, 2002). However, improved dating indicates that 
such periodicity is an artifact of ice core chronologies and not 
statistically significant (Ditlevsen et al., 2007), and inspec- 
tion of Fig. 27b reveals a broad range of timescales. Instead, 
the data imply a climate system that responds sensitively to 
even weak forcings and stochastic variability, both of which 
can spur amplifying feedbacks with a range of characteristic 
response times. 

Two conclusions are especially germane. First, subsur- 
face ocean warming is an effective mechanism for destabi- 
lizing ice shelves and thus the ice sheets buttressed by the 
ice shelves. Second, large rapid sea level rise can occur as a 
result of melting ice shelves. 

However, ice shelves probably were more extensive during 
glacial times. So are today’s ice sheets much more stable? 
The need to understand ice sheet vulnerability focuses atten- 
tion on end-Eemian events, when ice sheets were comparable 
in size to today’s ice sheets. 

4.2.4 End-Eemian climate and sea level change 

Termination II, ushering in the Eemian, was spurred by a late 
spring 60° N insolation anomaly peaking at +45 W m -2 at 

129.5 ky b2k (Fig. 27a), the largest anomaly in at least the 
past 425 ky (Fig. 3 of Hansen et al., 2007b). CO 2 and albedo 
forcings were mutually reinforcing. CO 2 began to rise before 
Antarctic <5 ls O, as deglaciation and warming began in the 
Northern Hemisphere. Most of the total CO 2 rise was pre- 
sumably from deep-ocean ventilation in the Southern Ocean, 
aided by meltwater that slowed the AMOC and thus helped 
to warm the Southern Ocean. 

The northern late spring insolation anomaly fell rapidly, 
becoming negative at 123.8 ky b2k (Fig. 27a), by which 
time summer insolation also began to fall (Fig. 26). North- 
ern Hemisphere ice sheets must have increased intermittently 
while Southern Hemisphere ice was still declining, consis- 
tent with minor, growing ice rafting events C27, C27a, C27b 
and C26 and a sea level minimum during 125-121 ky b2k 
(Sect. 4.1.1). High Eemian climate variability in the Antarc- 
tic (Pol et al., 2014) was likely a result of the see-saw relation 
with North Atlantic events. 
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C0 2 (Fig. 27c) remained at ~ 270 ppm for almost 15ky 
as the positive insolation anomaly on the Southern Ocean 
(Fig. 27a) kept the deep ocean ventilated. Sea level in the Red 
Sea analysis (Grant et ah, 2012) shown in Fig. 27f seems to 
be in decline through the Eemian, but that must be a combi- 
nation of dating and sea level error, as numerous sea level 
analyses cited in Sect. 4.1.1, our Supplement, and others 
(e.g., Chen et ah, 1991; Stirling et ah, 1998; Cutler et ah, 
2003) indicate high sea level throughout the Eemian and al- 
low a possible late-Eemian maximum. Chen et ah (1991), 
using a U-series dating with 2er uncertainty ±1.5 ky, found 
that the Eemian sea level high stand began between 132 and 
129 ky b2k, lasted for 12 ky, and was followed by rapid sea 
level fall. 

We assume that C26, the sharp cooling at 1 16.72 ky b2k 
in the NGRIP ice on the AICC2012 timescale, marks the end 
of fully interglacial Eemian conditions, described as 5e sensu 
stricto by Bauch and Erlenkeuser (2008). <i lx O in Antarctica 
was approaching a relative minimum (—46.7 %c at EDML; 
see Fig. S20 for details) and C0 2 was slowly declining at 
263 ppm. In the next 300 years 5 18 0 increased to —45.2 and 
C0 2 increased by 13 ppm with lag ~ 1500 years, which we 
interpret as see-saw warming of the Southern Ocean in re- 
sponse to the C26-induced AMOC slowdown and resulting 
increased SMOC ventilation of C0 2 . 

Freshwater causing the C26 AMOC shutdown could not 
have been Greenland surface melt. Greenland was already 
2000 years into a long cooling trend and the northern warm- 
season insolation anomaly was in the deepest minimum of 
the last 150 ky (Fig. 26a). Instead, C26 was one event in a se- 
ries, preceded by C27b and followed by C25, each a result of 
subsurface North Atlantic warming that melted ice shelves, 
causing ice sheets to discharge ice. Chapman and Shackle- 
ton (1999) did not find IRD from C26 in the mid- Atlantic, 
but Carlson et al. (2008) found a sharp increase in sediments 
near the southern tip of Greenland that they identified with 
C26. 

We suggest that the Southern Hemisphere was the source 
for brief late-Eemian sea level rise. The positive warm- 
season insolation anomaly on the Southern Ocean and 
AMOC slowdown due to C26 added to Southern Ocean heat, 
causing ice shelf melt, ice sheet discharge, and sea level rise. 
Rapid Antarctica ice loss would cool the Southern Ocean and 
increase sea ice cover, which may have left telltale evidence 
in ice cores. Indeed, Masson-Delmotte et al. (2011) suggest 
that abrupt changes of <5 ls O in the EDML and TALDICE ice 
cores (those most proximal to the coast) indicate a change 
in moisture origin, likely due to increased sea ice. Further 
analysis of Antarctic data for the late Eemian might help 
pinpoint the melting and help assess vulnerability of Antarc- 
tic ice sheets to ocean warming, but this likely will require 
higher-resolution models with more realistic sea ice distribu- 
tion and seasonal change than our present model produces. 

Terrestrial records in Northern Europe reveal rapid end- 
Eemian cooling. Sirocko et al. (2005) find cooling of 3 °C 
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in summer and 5-10 °C in winter in southern Germany, with 
annual layering in a dry Eifel maar lake revealing a 468-year 
period of aridity, dust storms, bushfires, and a decline of ther- 
mophilous trees. Similar cooling is found at other German 
sites and La Grande Pile in France (Kuhl and Litt, 2003). 
Authors in both cases interpret the changes as due to a south- 
ward shift of the polar front in the North Atlantic correspond- 
ing to C26. Cooling of this magnitude in northern Europe 
and increased aridity are found by Brayshaw et al. (2009) 
and Jackson et al. (2015) in simulations with high-resolution 
climate models forced by AMOC shutdown. 

While reiterating dating uncertainties, we note that the 
cool period with reduced NADW formation identified in re- 
cent high-resolution ocean core studies for Eirik Drift site 
MD03-2664 (Fig. 25) near Greenland (Irvali et al., 2012; 
Galaasen et al., 2014) at ~ 117 ky b2k has length similar to 
the 468-year cold stormy period found in a German lake core 
(Sirocko et al., 2005). 

The Eirik core data show a brief return to near-Eemian 
conditions and then a slow decline, similar to the oscillation 
in the NGRIP ice core at 1 16.72 ky b2k on the AICC2012 
timescale. 

The principal site of NADW formation may have moved 
from the GIN seas to just south of Greenland at end of the 
Eemian. Southward shift of NADW formation and the polar 
front is consistent with the sudden, large end-Eemian cooling 
in the North Atlantic and northern Europe, while cooling in 
southern European was delayed by a few millennia (Brauer et 
al., 2007). Thus end-Eemian midlatitude climate was charac- 
terized by an increased meridional temperature gradient, an 
important ingredient for strengthening storms. 

5 Modern data 

Observations help check our underlying assumption of non- 
linear meltwater growth and basic simulated climate effects. 
As these data are updated, and as more extensive observa- 
tions of the ocean and ice processes are obtained, a clearer 
picture should emerge over the next several years. 

5.1 Ice sheet mass loss and sea level rise 

The fundamental question we raise is whether ice sheet melt 
in response to rapid global warming will be nonlinear and 
better characterized by a doubling time for its rate of change 
or whether more linear processes dominate. Hansen (2005, 
2007) argued on heuristic grounds that ice sheet disintegra- 
tion is likely to be nonlinear if climate forcings continue to 
grow, and that sea level rise of several meters is possible on 
a timescale of the order of a century. Given current ice sheet 
melt rates, a 20-year doubling rate produces multi-meter sea 
level rise in a century, while 10- and 40-year doubling times 
require about 50 and 200 years, respectively. 
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Church et al. (2013) increased estimates of sea level rise 
compared to prior IPCC reports, but scenarios they dis- 
cuss are close to linear responses to the assumed rising cli- 
mate forcing. The most extreme climate forcing (RCP8.5, 
936 ppm CO 2 in 2100 and GHG forcing 8.5 Wm -2 ) is es- 
timated to produce 0.74 m sea level rise in 2100 relative to 
the 1986-2005 mean sea level, with the “likely” range of 
uncertainty at 0.52-0.98 m. Church et al. (2013) also dis- 
cuss semi-empirical estimates of sea level rise, which yield 
~ 0.7-1. 5 m for the RCP8.5 scenario, but express low confi- 
dence in the latter, thus giving preference to the model-based 
estimate of 0.52-0.98 m. We note that Sect. 4.4. 4.2 on ice 
sheet processes in the IPCC chapter on cryosphere observa- 
tions (Vaughan et al., 2013) contains valuable discussion of 
nonlinear ice sheet processes that could accelerate ice sheet 
mass loss but which are not fully included in current ice sheet 
models. 

Empirical analyses are needed if we doubt the realism of 
ice sheet models, but semi-empirical analyses lumping multi- 
ple processes together may yield a result that is too linear. Sea 
level rises as a warming ocean expands, as water storage on 
continents changes (e.g., in aquifers and behind dams), and 
as glaciers, small ice caps, and the Greenland and Antarctic 
ice sheets melt. We must isolate the ice sheet contribution, 
because only the ice sheets threaten multi-meter sea level 
rise. 

Hay et al. (2015) reanalyzed tide-gauge data for 1901- 
1990 including isostatic adjustment at each station, find- 
ing global sea level rise to be 1 .2 ± 0.2 mm year - 1 . Prior 
tide-gauge analyses of 1.6-1. 9 mm year -1 were inconsis- 
tent with estimates for each process, which did not add up 
to such a large value (Church et al., 2013). This estimate 
of 1.2 ±0.2 mm year -1 for 1900-1990 compares with esti- 
mated sea level rise of ~ 0.2 m in the prior two millennia or 
~ 0.1 mm year -1 (Kemp et al., 2011) and several estimates 
of~3 mm year - 1 for the satellite era (1993-present). Nerem 
et al. (2010) find sea level increase of 3.3 mm year -1 in the 
satellite era, while Watson et al. (2015), based in part on cal- 
ibration to tide-gauge data, suggest alternative rates of 2.9 or 
2.6 mm year -1 . 

Accepting the analyses of Hay et al. (2015) for 1901-1990 
and estimates of 2. 5-3. 5 mm year - 1 for the satellite era leads 
to a picture of a rising sea level rate (Fig. 29) that differs 
from the perception of near-linear sea level rise created by 
Fig. 13.3 in the IPCC report (Church et al., 2013). We do 
not argue for the details in Fig. 29 or suggest any change 
points for the rate of sea level rise, but the data do reveal a 
substantial increase in the rate of sea level rise. 

The majority of sea level rise in the 20th century was 
from the several processes other than Greenland and Antarc- 
tica mass loss (Church et al., 2013), so the timescale for ice 
sheet mass loss may differ from the timescale for past sea 
level change. A direct measure of ice sheet mass loss is ob- 
tained from satellite gravity measurements by Velicogna et 
al. (2014), who find Greenland’s mass loss in 2003-2013 
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Global mean sea level change 



Figure 29. Estimated sea level change (mm) since 1900. Data 
through 1992 are the tide-gauge record of Church and White (201 1) 
with the change rate multiplied by 0.78, so as to yield a mean 1901 — 
1990 change rate of 1.2 mm year -1 (Hay et al., 2015). The two esti- 
mates for the satellite era (1993-2015) are from Nerem et al. (2010, 
updated at http://sealevel.colorado.edu) and Watson et al. (2015). 

of280±58Gtyear -1 4 accelerating by 25.4 ± 1.2 Gtyear -2 , 
and Antarctic mass loss of 67 ± 44 Gt year - 1 accelerating by 
11 ±4Gtyear -2 . The Velicogna et al. (2014) data are up- 
dated in Fig. 30. The reduced mass loss rate of Greenland 
in 2013-2014 makes it difficult to infer an empirical growth 
rate for mass loss, as discussed below. 

Reliability of mass loss inferred from gravity data is 
supported by comparison to surface mass balance stud- 
ies (Fig. 30). Mass loss accelerations over 1992-2011 
obtained via the mass budget method (Rignot et al., 
2011) for Greenland (21.9 ± 1 Gtyear -2 ) and Antarctica 
(14.5 ± 2 Gt year -2 ) are similar to or larger than results from 
gravity analysis. A third approach, based on satellite radar 
altimetry, is consistent with the other two for mass loss from 
Greenland and West Antarctica (Shepherd et al., 2012), in- 
cluding the Amundsen Sea sector, which is the dominant 
contributor to Antarctic ice mass loss (Sutterley et al., 2014). 
Differences among techniques exist in East Antarctica, but 
mass changes there are small (Shepherd et al., 2012). 

Best-fit exponential doubling times for Greenland are 
4.8 years (MBM 1992-2010 data), 18.9 years (GRACE 
2003-2015 data) and 8.8 years (MBM + GRACE 1992- 
2015 data); in the latter case only GRACE data are used 
after 2002. The best fit to MBM ± GRACE is shown in 
Fig. 30. The equivalent results for Antarctica are 5.3, 3.2 and 


4 For comparison, our assumed freshwater injection of 
360 Gtyear -1 in 2011 with 10-year doubling yields an aver- 
age mass loss of 292 Gtyear -1 for 2003-2013. Further, Velicogna 
et al. (2014) find an ice mass loss of 74 ±7 Gtyear -1 from 
nearby Canadian glaciers and ice caps with acceleration of 
10 ± 2 Gtyear -2 , and there is an unknown freshwater input from 
melting ice shelves. Thus our assumed Northern Hemisphere 
meltwater was conservative. 


www.atmos-chem-phys.net/16/3761/2016/ 


Atmos. Chem. Phys., 16, 3761-3812, 2016 


3794 


J. Hansen et al.: Ice melt, sea level rise and superstorms 


(a) Greenland ice mass change rate (b) Antarctica ice mass change rate 



Figure 30. Greenland and Antarctic ice mass change. GRACE data are extension of Velicogna et al. (2014) gravity data. MBM (mass budget 
method) data are from Rignot et al. (2011). Red curves are gravity data for Greenland and Antarctica only; small Arctic ice caps and ice shelf 
melt add to freshwater input. 


5.0 years. Clearly the data records are too short to infer a 
doubling rate, let alone confirm that mass loss is exponen- 
tial. Thus we also show a 10-year doubling growth curve for 
Antarctic ice mass loss (Fig. 30b). The recent reduction of 
mass loss from Greenland illustrates how sensitive the empir- 
ical result is to record brevity, but the curves should become 
more informative over the next several years. 

Additional insight is provided by the regional breakdown 
of the mass change data as achieved in the Velicogna et 
al. (2014) analysis. The regional data suggest that the Antarc- 
tic situation may be more threatening than indicated by the 
continental mass loss rate. This net mass loss combines mass 
loss via ice streams with regions of net snow accumulation. 
Queen Maud Land, for example, is gaining 63 ± 6 Gt year 1 , 
accelerating by 15 ± 1 Gtyear 2 , but this mass gain may be 
temporary. Our simulations with increasing freshwater input 
indicate that circum-Antarctic cooling and sea ice increase 
eventually may limit precipitation reaching the continent, 
and recent SST and sea ice data have a tendency consistent 
with that expectation (Sect. 5.2). 

Amundsen Sea glaciers are a gateway to West Antarc- 
tic ice, which has potential for several meters of sea level. 
Mass loss of the Amundsen Sea sector was 1 1 6 ± 6 Gt year - 1 
in 2003-2013, growing 13 ± 2 Gtyear -2 (Velicogna et al., 
2014; Rignot et al., 2014; Sutterley et al., 2014). 

Totten Glacier in East Antarctica fronts the Aurora Sub- 
glacial Basin, which has the potential for ~ 6.7 m of sea level 
rise (Greenbaum et al., 2015). Williams et al. (2011) find 
that warm modified Circumpolar Deep Water is penetrating 
the continental shelf near Totten beneath colder surface lay- 
ers. Details of how warmer water reaches the ice shelf are 
uncertain (Khazendar et al., 2013), but, as in West Antarc- 
tica, the inland-sloping trough connecting the ocean with the 
main ice shelf cavity (Greenbaum et al., 2015) makes Tot- 
ten Glacier susceptible to unstable retreat (Goldberg et al., 
2009). Cook Glacier, further east in East Antarctica, also 
rests on a submarine inland-sloping bed and fronts ice equiv- 


alent to 3^1 m of sea level. The Velicogna et al. (2014) anal- 
ysis of gravity data for 2003-2013 finds the Totten sector of 
East Antarctica losing 17 ±4 Gtyear -1 , with the loss accel- 
erating by 4 ± 1 Gt year -2 , and the Victoria/Wilkes sector in- 
cluding Cook Glacier losing 16 ± 5 Gtyear -1 , with a small 
deceleration (2 ± 1 Gtyear -2 ). 

Greenland ice melt is subject to multiple feedbacks, some 
of which are largely absent on Antarctica, so it is not cer- 
tain whether Greenland ice is less or more vulnerable than 
Antarctic ice. On the one hand, some differences make the 
Greenland ice sheet seem less vulnerable. Greenland does 
not have as much unstable ice volume sitting behind retro- 
grade beds. Also, although surface cooling due to freshwater 
injection leads to subsurface ocean warming (Fig. 12), fresh- 
water injection may also reduce poleward transport of heat 
by the Atlantic Ocean if the AMOC slows down, and North 
Atlantic cooling may affect summer surface melt on Green- 
land. 

On the other hand, the Greenland ice sheet is subject to 
forcings and feedbacks that are less important on Antarc- 
tica. Greenland experiences extensive summer surface melt 
(Tedesco et al., 2011; Box et al., 2012), which makes sur- 
face albedo changes more important on Greenland. Green- 
land mass loss is thus more affected by snow darkening via 
dust, black carbon and biological substances including algae 
(Benning et al., 2014; Yasunari et al., 2015), which are in part 
an imposed climate forcing but in some cases also a substan- 
tial amplifying feedback. Soot from forest fires occurs natu- 
rally, but the magnitude of fire events is increasing (Flanni- 
gan et al., 2013; Jolly et al., 2015) and may have contributed 
to widespread Greenland melt events in recent years (Keegan 
et al., 2014). Pigmented algae can substantially reduce spring 
and summer ice albedo and may be an important feedback 
in a warming world (Benning et al., 2014). Other amplify- 
ing feedbacks for Greenland include the ice surface eleva- 
tion feedback, cryo-hydrologic warming in which percolat- 
ing water alters thermal regime and weakens the ice sheet on 
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decadal timescales (Colgan et al., 2015) and ocean-mediated 
melting of ice shelves and glacier fronts (Rignot et al., 2010). 
Increasing ice sheet surface melt and increasing ice stream 
mass discharge are both contributing to the observed growing 
mass loss rate of the Greenland ice sheet, as discussed in our 
response AC7962 on the ACPD website. Such mutually rein- 
forcing processes provide an expectation of nonlinear mass 
loss increase if the climate forcing continues to increase. 

Interpretation of Greenland mass loss is made difficult by 
its high variability. Large 2010-2012 mass loss was related 
to unusual summer high pressure over Greenland (Fettweis 
et al., 2013; Bellflamme et al., 2015), which produced a per- 
sistent “atmospheric river” of warm air of continental origin 
(Neff et al., 2014). However, weather patterns were much less 
favorable for surface melt in 2013 and 2014, and Greenland 
mass loss was much reduced (Fig. 30a). 

We conclude that empirical data are too brief to imply a 
characteristic time for ice sheet mass loss or to confirm our 
hypothesis that continued high fossil fuel emissions leading 
to CO2 ~ 600-900 ppm will cause exponential ice mass loss 
up to several meters of sea level. The empirical data are con- 
sistent with a doubling time of the order of a decade, but they 
cannot exclude slower responses. Our expectation of nonlin- 
ear behavior is based in part on recognition of how multiple 
amplifying feedbacks feed upon each other (Hansen et al., 
1984; Pollard et al., 2015) and thus can result in large rapid 
change. 

5.2 Sea surface temperature and sea ice 

The fundamental difference between climate forecasts of 
our model and CMIP simulations employed in IPCC assess- 
ments should appear in sea surface temperature (SST), which 
is well monitored. The Southern Ocean warms steadily in 
CMIP5 models (Fig. 12.11 of Collins et al., 2013) that have 
little or no freshwater injection. In contrast, ice melt causes 
Southern Ocean cooling in our model, especially in the West- 
ern Hemisphere (Fig. 16). The model’s cooling is largest in 
the Western Hemisphere because the specified freshwater in- 
jection (Fig. 14), based on data of Rignot et al. (2013) and 
Depoorter et al. (2013), is largest there. The cooling pat- 
tern is very strong by 2055-2060 (Fig. 16), when freshwa- 
ter injection reaches 3.8 Sv years on the North Atlantic and 
7.6 Sv years on the Southern Ocean, amounts that yield 1 m 
global sea level rise with one-third from Northern Hemi- 
sphere ice. SST observations already show a cooling trend 
in the Southern Ocean off West Antarctica (Fig. 31) and are 
suggestive that the real world may be more sensitive than 
the model, but additional years of data are needed to confirm 
that. 

Our model also differs from models in the predicted sense 
of Southern Hemisphere sea ice change. Freshwater effects 
dominate over direct effects of GHGs in our model, and thus 
sea ice cover grows. Thompson et al. (201 1) suggest that O3 
depletion may account for observed Antarctic sea ice growth. 
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but Sigmond and Fyfe (2014) found that all CMIP5 models 
yield decreasing sea ice in response to observed changes of 
O3 and other GHGs. Ferreira et al. (2015) show that O3 de- 
pletion yields a short timescale sea ice increase that is soon 
overtaken by warming and sea ice decrease with realistic 
GHG forcing. We suggest that these models are missing the 
dominant driver of change on the Southern Ocean: freshwa- 
ter input. 

Our modeled SMOC has begun to slow already (Fig. 32a), 
consistent with tracer observations in the Weddell Sea by 
Huhn et al. (2013), which reveal a 15-21 % reduction in the 
ventilation of Weddell Sea Bottom Water and Weddell Sea 
Deep Water in 1984—2008. Delayed growth of sea ice in the 
model (Fig. 32b) may be in part related to the model’s muted 
vertical stratification, as we will discuss, and the model’s 
general difficulty in producing Southern Hemisphere sea ice. 

We infer that observed cooling in the western part of the 
Southern Ocean, growing Southern Ocean sea ice, and slow- 
down of at least the Weddell Sea component of SMOC are 
early responses to increasing freshwater injection. Although 
observed sea ice increase is smaller in 2015 than in the previ- 
ous few years (data are updated daily at http://nsidc.org/data/ 
seaice_index), but Hansen and Sato (2016) note a negative 
correlation between sea ice area and El Ninos, so we expect 
that sea ice growth may resume after the present strong El 
Nino fades. 

Let us compare North Atlantic and Southern Ocean re- 
sponses to freshwater forcing. The modeled AMOC response 
does not become significant until ~ 2040 (Fig. 33). Even dis- 
counting the decade lead in Southern Hemisphere forcing 
(720 Gt year -1 in 2011, double that in the North Atlantic), 
SMOC still responds quicker, albeit gradually, to freshwater 
forcing (Fig. 32). 

Observations suggest that the real world may be respond- 
ing more quickly than the model to freshwater forcing in the 
North Atlantic. Rahmstorf et al. (2015) develop an AMOC 
index based on SST in the “global warming hole” south- 
east of Greenland (Drijfhout et al., 2012), and they use the 
AMOC index to conclude that an AMOC slowdown un- 
precedented in the prior 1000 years occurred in the late 20th 
century. That slowdown seems to have been a response to 
the “Great Salinity Anomaly”, which is thought to have re- 
sulted from an estimated ~ 2000 km 3 anomalous sea ice ex- 
port from the Arctic (Dickson et al., 1988). Although the 
AMOC partially recovered in the early 21st century, fur- 
ther slowdown has returned in the past several years, judging 
from a measurement array (Robson et al., 2014) as well as 
from the AMOC index (Rahmstorf et al., 2015). The recent 
AMOC slowdown could be related to ice melt, as Greenland 
(Fig. 30) and neighboring ice caps contributed more than 
1000 km 3 meltwater in 2011-2012. The model (Fig. 33), in 
contrast, does not reach substantial AMOC reduction until 
~ 2040, when the annual freshwater injection into the North 
Atlantic is ~ 7500 km 3 year -1 and the cumulative injection 
is ~ 1.2 Sv years. 
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(a) SST change in past 25 years ( °C) 


1990-2015 


0.30 



(b) West Antarctic SST with 1979-2000 base period 



Figure 31. (a) Observed 1990-2015 SST change based on local linear trends and (b) SST anomaly relative to 1979-2000 for area south of 
56° S between the dateline and 50° W. Base period excludes data prior to 1979 to avoid use of Southern Ocean climatology that artificially 
reduces variability (Huang et al., 2015). 


(a) SMOC at 72° S, 2. 1 km depth (b) Southern hemispheric mean sea ice area anomaly 



Figure 32. (a) Global meridional overturning circulation (Sv) at 72° S. Freshwater injection near Antarctica is 720Gtyear~ l in 2011, 
increasing with 10-year doubling time, and half as much around Greenland. SMOC diagnostic includes only the mean (Eulerian) term, 
(b) Annual-mean Southern Hemisphere sea ice area anomaly (10 6 km 2 ) in five runs (relative to 1979-2000). Observations include 2015. 


A useful calibration of AMOC sensitivity to freshwater 
forcing is provided by the 8.2 ky b2k glacial Lake Agassiz 
freshwater outburst (Kleiven et al., 2008) that occurred with 
the demise of the Hudson Bay ice dome. Freshwater injected 
into the North Atlantic was ~ 2.5-5 Sv years (Clarke et al., 
2004). Proxy temperature records (see LeGrande et al., 2006) 
suggest that real-world cooling reached about 6 °C on Green- 
land, 3^f°C in the east Norwegian Sea and 1.7 °C in Ger- 
many. The duration of the 8.2 ky b2k event was 160 years 
(Rasmussen et al., 2014). The LeGrande et al. (2006) model, 
which has the same atmospheric model as our present model 
but does not include the basic improvements in the ocean de- 
scribed in Sect. 3.2, produced results not inconsistent with 
real-world changes, but the modeled temperature response 
seemed to be on the low side (Fig. 1 of LeGrande et al., 
2006). The mean decrease in the modeled AMOC was 40 %, 
with AMOC recovering within 20-30 years but secondary 
and tertiary slowdowns in some of the model runs extending 
as long as in the observed 8.2 ky b2k event (160 years). Al- 
though this model response is within the range suggested by 
paleo-data, it is on the weak side. 

This model check based on the 8.2 ky b2k event does 
not prove that the model has correct sensitivity for today’s 
weaker forcing. We suspect that the model is less sensitive 
than the real world because the model has difficulty main- 
taining vertical stratification, which could result from coarse 


vertical resolution, excess parameterized small-scale mixing, 
or numerical noise. Excessive mixing could also explain a 
too-long climate response time, as discussed in connection 
with Fig. 4. Hansen et al. (201 1) showed that surface temper- 
ature is probably too sluggish in response to a climate forc- 
ing, not only in the GISS model but also in several other mod- 
els. Hofmann and Rahmstorf (2009) suggest another reason 
for models being too insensitive to freshwater forcing: a bias 
in ocean model development spurred by desire for a stable 
AMOC. Below we suggest studies that are needed to investi- 
gate the model sensitivity issue. 

5.3 Southern Ocean internal processes 

Although the ocean surface is observed in detail on a daily 
basis, our main interest is in implications for long-term 
processes in the ocean below. Paleoclimate data discussed 
in Sect. 4 reveal that the Southern Ocean, as a gateway 
to the global deep ocean, exerts a powerful control over 
glacial/interglacial climate. 

The Southern Ocean has significant control on release of 
ocean heat to space. In an extreme case, poly ny as form in 
the dead of Antarctic winter, as upwelling warm water melts 
the sea ice and raises the air temperature by tens of degrees, 
increasing thermal radiation to space, thus serving as a valve 
that releases ocean heat. Today, as surface meltwater stabi- 
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(a) AMOC at 28° N 


(b) SST in "global warming hole" region 



Figure 33. (a) AMOC (Sv) at 28° N in simulations with the forcings of Sect. 4.2 (i.e., including freshwater injection of 720 Gt year 1 in 
2011 around Antarctica, increasing with a 10-year doubling time, and half that amount around Greenland), (b) SST (°C) in the North Atlantic 
region (44-60° N, 10-50° W). 


lizes the vertical water column, that valve is being partially 
closed. De Lavergne et al. (2014) relate the absence of large 
open-ocean polynyas in recent decades to surface freshening. 
Release of heat to the atmosphere and space, which occurs 
without the need for large open-ocean polynyas, is slowed 
by increasing sea ice cover in response to increasing ice shelf 
melt (Bintanja et al., 2013). 

Internal Southern Ocean effects of ocean surface fresh- 
ening and cooling seem to be well underway. Schmidtko et 
al. (2014) and Roemmich et al. (2015) document changes in 
the Southern Ocean in recent decades, especially warming 
of Circumpolar Deep Water (CDW), which they and others 
(Jacobs et al., 2011; Rignot et al., 2013) note is the likely 
cause of increased ice shelf melt. Observations of ocean sur- 
face freshening and freshening of the water column (Rintoul, 
2007; Jacobs and Giulivi, 2010) and deep-ocean warming 
(Johnson et al., 2007; Purkey and Johnson, 2013) leave lit- 
tle doubt that these processes are occurring. 

6 Summary and implications 

Via a combination of climate modeling, paleoclimate anal- 
yses, and modern observations we have identified climate 
feedback processes that help explain paleoclimate change 
and may be of critical importance in projections of human- 
made climate change. Here we summarize our interpretation 
of these processes, their effect on past climate change, and 
their impact on climate projections. We then discuss key ob- 
servations and modeling studies needed to assess the validity 
of these interpretations. We argue that these feedback pro- 
cesses may be understated in our model, and perhaps other 
models, because of an excessively diffusive ocean model. 
Thus, there is urgency to obtain a better understanding of 
these processes and models. 


6.1 Ocean stratification and ocean warming 

Global ocean circulation (Fig. 16) is altered by the effect of 
low-density freshwater from melting of Greenland or Antarc- 
tic ice sheets. While the effects of shutdown of NADW have 
been the subject of intensive research for a quarter of a cen- 
tury, we present evidence that models have understated the 
threat and imminence of slowdown and shutdown AMOC 
and SMOC. Below we suggest modeling and observations 
that would help verify the reality of stratification effects on 
polar oceans and improve assessment of likely near-term and 
far-term impacts. 

Almost counter-intuitively, regional cooling from ice melt 
produces an amplifying feedback that accelerates ice melt by 
placing a lid on the polar ocean that limits heat loss to the 
atmosphere and space, warming the ocean at the depth of ice 
shelves. The regional surface cooling increases Earth’s en- 
ergy imbalance, thus pumping into the ocean energy required 
for ice melt. 5 

6.2 Southern Ocean, CO 2 control knob, and ice sheet 
timescale 

Our climate simulations and analysis of paleoclimate oscilla- 
tions indicate that the Southern Ocean has the leading role in 
global climate change, with the North Atlantic a supporting 
actor. The Southern Ocean dominates by controlling ventila- 
tion of the deep-ocean CO 2 reservoir. 

CO 2 is the control knob that regulates global tempera- 
ture. On short timescales, i.e., fixed surface climate, CO 2 
sets atmospheric temperature because CO 2 is stable; thus, the 
ephemeral radiative constituents, H 2 O and clouds, adjust to 
CO 2 amount (Lacis et al., 2010, 2013). 


5 Planetary energy imbalance induced by meltwater cooling 
helps provide the energy required by ice heat of fusion. Ice melt 
to raise sea level 1 m requires a 10-year Earth energy imbalance 
0.9 Wm“ 2 (Table SI; Hansen et al., 2005b). 
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On millennial timescales both CO 2 and surface albedo (de- 
termined by ice and snow cover) are variable and contribute 
about equally to global temperature change (Hansen et al., 
2008). However, here too CO 2 is the more stable constituent 
with timescale for change ~ 10 3 years, while surface albedo 
is more ephemeral judging from the difficulty of finding any 
lag of more than the order of 10 2 years between sea level and 
polar temperature (Grant et al., 2012). 

Here we must clarify that ice and snow cover are both 
a consequence of global temperature change, generally re- 
sponding to the CO 2 control knob, but also a mechanism for 
global climate change. Specifically, regional or hemispheric 
snow and ice respond to seasonal insolation anomalies (as 
well as to CO 2 amount), thus affecting hemispheric and 
global climate, but to achieve large global change the albedo- 
driven climate change needs to affect the CO 2 amount. 

We also note that Southern Ocean ventilation is not the 
only mechanism affecting airborne CO 2 amount. Terrestrial 
sources, dust fertilization of the ocean, and other factors play 
roles, but deep-ocean ventilation seems to be the dominant 
mechanism on glacial-interglacial timescales. 

The most important practical implication of this “control 
knob” analysis is realization that the timescale for ice sheet 
change in Earth’s natural history has been set by CO 2 , not by 
ice physics. With the rapid large increase in CO 2 expected 
this century, we have no assurance that large ice sheet re- 
sponse will not occur on the century timescale or even faster. 

6.3 Heinrich and Dansgaard-Oeschger events 

Heinrich and Dansgaard-Oeschger events demonstrate the 
key role of subsurface ocean warming in melting ice shelves 
and destabilizing ice sheets, and they show that melting ice 
shelves can result in large rapid sea level rise. A cold lens of 
fresh meltwater on the ocean surface may make surface cli- 
mate uncomfortable for humans, but it abets the provision of 
warmth at depths needed to accelerate ice melt. 

6.4 End-Eemian climate events 

We presented evidence for a rapid sea level rise of several 
meters late in the Eemian, as well as evidence of extreme 
storms in the Bahamas and Bermuda that must have occurred 
when sea level was near its maximum. This evidence is con- 
sistent with the fact that the North Atlantic was cooling in the 
late Eemian, while the tropics were unusually warm, the lat- 
ter being consistent with the small obliquity of Earth’s spin 
axis at that time. 

Giant boulders of mid-Pleistocene limestone placed atop 
an Eemian substrate in North Eleuthera, which must have 
been deposited by waves, are emblematic of stormy end- 
Eemian conditions. Although others have suggested the boul- 
ders may have been emplaced by a tsunami, we argue that 
the most straightforward interpretation of all evidence favors 
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storm emplacement. In any case, there is abundant evidence 
for strong late-Eemian storminess and high sea level. 

A late-Eemian shutdown of the AMOC would have caused 
the most extreme North Atlantic temperature gradients. 
AMOC shutdown in turn would have added to Southern 
Ocean warmth, which may have been a major factor in the 
Antarctic ice sheet collapse that is required to account for 
the several meters of rapid late-Eemian sea level rise. 

Confirmation of the exact sequence of late-Eemian events 
does not require absolute dating, but it probably requires 
finding markers that allow accurate correlation of high- 
resolution ocean cores with ice cores, as has proved possible 
for correlating Antarctic and Greenland ice cores. Such ac- 
curate relative dating would make it easier to interpret the 
significance of abrupt changes in two Antarctic ice cores 
at about end-Eemian time (Masson-Delmotte et al., 2011), 
which may indicate rapid large change in Antarctic sea ice 
cover. 

Understanding end-Eemian storminess is important in part 
because the combination of strong storms with sea level rise 
poses a special threat. However, sea level rise itself is the 
single greatest global concern, and it is now broadly ac- 
cepted that late-Eemian sea level reached +6-9 m, implicat- 
ing a substantial contribution from Antarctica, at a time when 
Earth was little warmer than today (Dutton et al., 2015; Sup- 
plement to our present paper). 

6.5 The Anthropocene 

The Anthropocene (Crutzen and Stoermer, 2000), the era in 
which humans have contributed to global climate change, 
is usually assumed to have begun in the past few centuries. 
Ruddiman (2003) suggests that it began earlier, as deforesta- 
tion began to affect CO 2 about 8000 years ago. Southern 
Ocean feedbacks considered in our present paper are relevant 
to that discussion. 

Ruddiman (2003) assumed that 40 ppm of human-made 
CO 2 was needed to explain a 20 ppm CO 2 increase in the 
Holocene (Fig. 27c), because CO 2 decreased by ~ 20 ppm, 
on average, during several prior interglacials. Such a large 
human source should have left an imprint on <5 13 C02 that 
is not observed in ice core CO 2 (Elsig et al., 2009). Ruddi- 
man (2013) suggests that 13 C was taken up in peat formation, 
but the required peat formation would be large and no persua- 
sive evidence has been presented to support such a dominant 
role for peat in the glacial carbon cycle. 

We suggest that Ruddiman’s hypothesis may be right, but 
the required human-made carbon source is much smaller than 
he assumed. Decline of CO 2 in interglacial periods is a cli- 
mate feedback, a result of declining Southern Ocean tem- 
perature, which slows the ventilation of the deep ocean and 
exhalation of deep-ocean CO 2 . Human-made CO 2 forcing 
needed to avoid Antarctic cooling and atmospheric CO 2 de- 
cline is only the amount needed to counteract the weak natu- 
ral forcing trend, not the larger feedback-driven CO 2 declines 
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in prior interglacials, because the feedback does not occur if 
the natural forcings are counteracted. 

The warm-season insolation anomaly on the Southern 
Ocean was positive and growing 8 ky ago (Fig. 27a). Thus the 
human-made CO 2 contribution required to make the South- 
ern Ocean a CO 2 source sufficient to yield the observed 
CO 2 growth (Fig. 27c) is unlikely to have been larger than 
~ 10 ppm, but quantification requires carbon cycle modeling 
beyond present capabilities. 

However, the modest requirement on the human CO 2 
source and the low 5 13 C content of deep-ocean CO 2 make 
the Ruddiman hypothesis more plausible and likely. 

6.6 The “Hyper-Anthropocene” 

A fundamentally different climate phase, a “Hyper- 
Anthropocene”, began in the latter half of the 18th century 
as improvements of the steam engine ushered in the indus- 
trial revolution (Hills, 1993) and exponential growth of fos- 
sil fuel use. Human-made climate forcings now overwhelm 
natural forcings. CO 2 , at 400 ppm in 2015, is off the scale in 
Fig. 27c. CO 2 climate forcing is a reasonable approximation 
of the net human forcing, because forcing by other GHGs 
tends to offset negative human forcings, mainly aerosols 
(Myhre et al., 2013). Most of the CO 2 growth occurred in the 
past several decades, and three-quarters of the ~ 1 °C global 
warming since 1850 (update of Hansen et al., 2010, available 
at http://www.columbia.edu/~mhsll9/Temperature/) has oc- 
curred since 1975. Climate response to this CO 2 level, so far, 
is only partial. 

Our analysis paints a very different picture than 
IPCC (2013) for continuation of this Hyper-Anthropocene 
phase, if GHG emissions continue to grow. In that case, we 
conclude that multi-meter sea level rise would become prac- 
tically unavoidable, probably within 50-150 years. Full shut- 
down of the North Atlantic Overturning Circulation would be 
likely within the next several decades in such a climate forc- 
ing scenario. Social disruption and economic consequences 
of such large sea level rise, and the attendant increases in 
storms and climate extremes, could be devastating. It is not 
difficult to imagine that conflicts arising from forced mi- 
grations and economic collapse might make the planet un- 
governable, threatening the fabric of civilization. 

Our study, albeit with a coarse-resolution model and sim- 
plifying assumptions, raises fundamental questions that point 
toward specific modeling and measurement needs. 

6.7 Modeling priorities 

Predictions from our modeling are shown vividly in Fig. 16, 
which shows simulated climate four decades in the future. 
However, we concluded that the basic features there are al- 
ready beginning to evolve in the real world, that our model 
underestimates sensitivity to freshwater forcing and the strat- 
ification feedback, and that the surface climate effects are 
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likely to emerge sooner than models suggest, if GHG climate 
forcing continues to grow. 

This interpretation arises from evidence of excessive 
small-scale mixing in our ocean model and some other mod- 
els, which reduces the stratification feedback effect of fresh- 
water injection. Our climate model, with ~ 3 °C equilibrium 
sensitivity for 2 x CO 2 , achieves only ~ 60 % of its equilib- 
rium response in 100 years (Fig. 4). Hansen et al. (201 1) con- 
clude that such a slow response is inconsistent with Earth’s 
measured energy imbalance; if the ocean were that diffusive 
it would be soaking up heat faster than the measured plan- 
etary energy imbalance ~0.6Wm -2 (Hansen et al., 2011; 
von Schuckmann et al., 2016). Hansen (2008) found the re- 
sponse time of climate models of three other modeling cen- 
ters to be as slow as or slower than the GISS model, im- 
plying that the oceans in those models were also too dif- 
fusive and thus their climate response times too long. The 
climate response time is fundamental to interpretation of cli- 
mate change and the impact of excessive small-scale mixing, 
if such exists, is so important that we suggest that all mod- 
eling groups participating in future CMIP studies should be 
asked to calculate and report their climate response function, 
R (Fig. 4). An added merit of that information is the fact 
that R permits easy calculation of the global temperature re- 
sponse for any climate forcing (Eq. 1). 

It may be possible to quickly resolve or at least clarify this 
modeling issue. A fundamental difficulty with ocean mod- 
eling is that the scale of eddies and jet-like flows is much 
smaller than comparable features in the atmosphere, which 
is the reason for the Gent and McWilliams (1990) parame- 
terization of eddy mixing in coarse-resolution models. How- 
ever, computers at large modeling centers today allow simu- 
lations with ocean resolution as fine as ~ 0.1°, which can re- 
solve eddies and minimize need for parameterizations. Win- 
ton et al. (2014) used a GFDL model (one of the models 
Hansen, 2008, found to have R similar to that of our model) 
with 0.1° ocean resolution for a 1 %year -1 increasing CO 2 
experiment, finding an ~ 25 % increase in transient global 
warming, which is about the increment needed to increase R 
(100 years) to ~0.75, consistent with Earth’s measured en- 
ergy imbalance (Hansen et al., 2011). The increased surface 
response implies that small-scale mixing that limits stratifi- 
cation is reduced. Saba et al. (2016) show that this model 
with 0.1° ocean resolution yields 3-4 °C warming along 
the United States East Coast at doubled CO 2 and cooling 
(~ — 1 °C) southeast of Greenland, both temperature changes 
a result of AMOC slowdown that reduces poleward transport 
of heat. 

The model results are striking because similar temperature 
patterns seem to be emerging in observations (Figs. 31, S24). 
Annual and decadal variability limit interpretation, but given 
the AMOC sensitivity revealed in paleoclimate data, we infer 
that stratification effects are beginning to appear in the North 
Atlantic due to the combination of ice melt and GHG forc- 
ing. Eddy-resolving ocean models are just beginning to be 
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employed and analyzed (Bryan et al., 2014), but there needs 
to be an added focus in CMIP runs to include freshwater 
from ice melt. CMIP5 simulations led to IPCC estimates of 
AMOC weakening in 2100 (Collins et al., 2013) of only 1 1 % 
for the weakest forcing scenario and 34 % for the strongest 
forcing (CO 2 = 936 ppm), but the CMIP5 runs do not in- 
clude ice melt. This moderate change on a century timescale 
may be a figment of (1) excluding ice melt and (2) under- 
stated stratification, as can be checked with improved high- 
resolution models that include realistic meltwater injection. 
Reliable projections of AMOC and North Atlantic climate 
will not flow simply from new high-resolution model runs, 
as Winton et al. (2014) note that ocean models have other 
tuning parameters that can sensitively affect AMOC stability 
(Hofmann and Rahmstorf, 2009), which is reason for a broad 
comparative study with the full set of CMIP models. 

High resolution ocean models are also needed to realisti- 
cally portray deepwater formation around Antarctica, pen- 
etration of warm waters into ice shelf environments, and, 
eventually, ocean-ice sheet feedbacks. More detailed mod- 
els should also include the cooling effect of ice phase change 
(heat of fusion) more precisely, perhaps including iceberg 
tracks. However, there is merit in also having a coarser- 
resolution version of major models with basically the same 
model physics. Coarser resolution allows long simulations, 
facilitating analysis of the equilibrium response, paleocli- 
mate studies, and extensive testing of physical processes. 

6.8 Measurement priorities 

A principal issue is whether ice melt will increase expo- 
nentially, as we hypothesize if GHGs continue to grow 
rapidly. Continuous gravity measurements, coupled with sur- 
face mass balance and physical process studies on the ice 
sheets, are needed to obtain and understand regional ice mass 
loss on both Greenland and Antarctica. Ocean-ice shelf in- 
teractions need to be monitored, especially in Antarctica, but 
some Greenland ice is also vulnerable to thermal forcing by 
a warming ocean via submarine glacial valleys (Morlighem 
et al., 2014; Khan et al., 2014). 

Summer weather variability makes mass loss in the Green- 
land melt season highly variable, but continued warming of 
North American continental air masses likely will spur mul- 
tiple amplifying feedbacks. These feedbacks need to be mon- 
itored and quantified because their combination can lead to 
rapid meltwater increase. In addition to feedbacks discussed 
in Sect. 5.1, Machguth et al. (2016) note that meltwater in- 
jection to the ocean will increase as surface melt and refreeze 
limits the ability of firn to store meltwater. Meltwater in the 
past several years is already of the magnitude of the “Great 
Salinity Anomaly” (Dickson et al., 1988) that Rahmstorf et 
al. (2015) conclude produced significant AMOC slowdown 
in the late 20th century. 

Continued global measurements of SST from satellites, 
calibrated with buoy and ship data (Huang et al., 2015) will 
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reveal whether coolings in the Southern Ocean and south- 
east of Greenland are growing. Internal ocean temperature, 
salinity and current measurements by the ARGO float pro- 
gram (von Schuckmann et al., 2016), including planned ex- 
tensions into the deep ocean and under sea ice, are crucial 
for several reasons. ARGO provides global measurements 
of ocean quantities that are needed to understand observed 
surface changes. If climate models are less sensitive to sur- 
face forcings than the real world, as we have concluded, the 
ARGO data will help analyze the reasons for model short- 
comings. In addition, ARGO measurements of the rate of 
ocean heat content change are the essential data for accu- 
rate determination of Earth’s energy imbalance, which de- 
termines the amount of global warming that is still “in the 
pipeline” and the changes of atmospheric composition that 
would be needed to restore energy balance, the fundamental 
requirement for approximately stabilizing climate. 

6.9 Practical implications 

The United Nations Framework Convention on Climate 
Change (UNFCCC, 1992) states the following: 

The ultimate objective of this Convention and any 
related legal instruments that the Conference of the 
Parties may adopt is to achieve, in accordance with 
the relevant provisions of the Convention, stabi- 
lization of greenhouse gas concentrations in the at- 
mosphere at a level that would prevent dangerous 
anthropogenic interference with the climate sys- 
tem. Such a level should be achieved within a time 
frame sufficient to allow ecosystems to adapt nat- 
urally to climate change, to ensure that food pro- 
duction is not threatened and to enable economic 
development to proceed in a sustainable manner. 

“Dangerous” is not further defined by the UNFCCC. Our 
present paper has several implications with regard to the con- 
cerns that the UNFCCC is meant to address. 

First, our conclusions suggest that a target of limiting 
global warming to 2 °C, which has sometimes been dis- 
cussed, does not provide safety. We cannot be certain that 
multi-meter sea level rise will occur if we allow global warm- 
ing of 2 °C. However, we know the warming would remain 
present for many centuries, if we allow it to occur (Solomon 
et al., 2010), a period exceeding the ice sheet response time 
implied by paleoclimate data. Sea level reached +6-9 m in 
the Eemian, a time that we have concluded was probably 
no more than a few tenths of a degree warmer than today. 
We observe accelerating mass losses from the Greenland and 
Antarctic ice sheets, and we have identified amplifying feed- 
backs that will increase the rates of change. We also observe 
changes occurring in the North Atlantic and Southern oceans, 
changes that we can attribute to ongoing warming and ice 
melt, which imply that this human-driven climate change 
seems poised to affect these most powerful overturning ocean 
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circulation systems, systems that we know have had huge ef- 
fects on the planetary environment in the past. We conclude 
that, in the common meaning of the word danger, 2 °C global 
warming is dangerous. 

Second, our study suggests that global surface air temper- 
ature, although an important diagnostic, is a flawed metric 
of planetary “health”, because faster ice melt has a cooling 
effect for a substantial period. Earth’s energy imbalance is 
in some sense a more fundamental climate diagnostic. Sta- 
bilizing climate, to first order, requires restoring planetary 
energy balance. The UNFCCC never mentions temperature 
- instead it mentions stabilization of greenhouse gas concen- 
trations at a level to avoid danger. It has been shown that the 
dominant climate forcing, CO 2 , must be reduced to no more 
than 350 ppm to restore planetary energy balance (Hansen 
et al., 2008) and keep climate near the Holocene level, if 
other forcings remain unchanged. Rapid phasedown of fossil 
fuel emissions is the crucial need, because of the millennial 
timescale of this carbon in the climate system. Improved un- 
derstanding of the carbon cycle is needed to determine the 
most effective complementary actions. It may be feasible to 
restore planetary energy balance via improved agricultural 
and forestry practices and other actions to draw down at- 
mospheric CO 2 amount, if fossil fuel emissions are rapidly 
phased out. 

Third, not only do we see evidence of changes beginning 
to happen in the climate system, as discussed above, but we 
have also associated these changes with amplifying feedback 
processes. We understand that in a system that is out of equi- 
librium, a system in which the equilibrium is difficult to re- 
store rapidly, a system in which major components such as 
the ocean and ice sheets have great inertia but are begin- 
ning to change, the existence of such amplifying feedbacks 
presents a situation of great concern. There is a possibility, a 
real danger, that we will hand young people and future gener- 
ations a climate system that is practically out of their control. 

We conclude that the message our climate science deliv- 
ers to society, policymakers, and the public alike is this: we 
have a global emergency. Fossil fuel CO 2 emissions should 
be reduced as rapidly as practical. 

The Supplement related to this article is available online 
at doi:10.5194/acp-16-3761-2016-supplement. 
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